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Abstract. We consider the spectral problem for the two-dimensional Schrodinger operator for a charged par- 
ticle in strong uniform magnetic and periodic electric fields. The related classical problem is analyzed first by 
means of the Krylov-Bogoljubov-Alfven and Neishtadt averaging methods. It allows us to show "almost inte- 
' grability" of the the original two-dimensional classical Hamilton system, and to reduce it to a one-dimensional 

one on the phase space which is a two-dimensional torus. Using the topological methods for integrable Hamil- 
tonian system and elementary facts from the Morse theory, we give a general classification of the classical 
motion. According this classification the classical motion is separated into different regimes with different 
topological characteristics (like rotation numbers and Maslov indices). Using these regimes, the semiclassi- 
• cal approximation, the Bohr-Sommerfeld rule and the correspondence principle, we give a general asymptotic 

description of the (band) spectrum of the original Schrodinger operator and, in particular, estimation for the 
number of subbands in each Landau band. From this point of view the regimes, are the classical preimages of 
^ ' "spectral series" of the Schrodinger operator. We also discuss the relationship between this spectrum and the 

04 ' spectrum of one-dimensional difference operators. 
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Both classical and quantum problems describing the motions of particles under the influence of a uniform 
'. magnetic and a periodic electric fields have very curious properties even in two dimensions. This has caused 

a large number of publications; which we mention here only some of them which are most relevant to our 



o 



^ . considerations [6-8,11,14,32,37,38,46,49-51,53,57,60,62,64,73-77,83,90,94-98]. 

^1 If the magnetic field is strong enough, then a large parameter appears in both classical and quantum 

^ ■ mechanics. Hence it is possible to use the averaging methods [2,3, 16-18,63,70,71,78] and semiclassical 

. approximation [9, 28, 47,48, 67-69]. Even though this circle of problems is well studied, we propose here 

^ I some apparently new formulas and interpretations. 

> 

X. 

^ , 1.1. The two-dimensional magnetic Schrodinger operator in a periodic electric field. Assumptions 
and parameters. We want to describe certain asymptotic spectral properties of the Schrodinger operator 



1. Formulation of the problem and brief description of the results 



(1.1) 



in L^(M^), which is essentially self-adjoint on (7q°(M^), as /i, e 0. We assume that the potential X2) 
is real analytic in and periodic with respect to the lattice F generated by two linearly independent vectors 
ai = (an, 012) = (2vr, 0), 02 = (021, 022), i- e. we have v{x + ai) = v{x + 02) = v{x). 

Such a problem arises in the following physical situation. Consider the motion of a particle with charge 
— e and mass m in the plane M.y in a uniform magnetic and periodic electric field. If the magnetic field is 
perpendicular to the y-plane and has strength B > 0, then this motion is described (in the Landau gauge) 
by the operator 



where V is the potential of the electric field and c, h are physical constants. Let the potential V is periodic 
with respect to the lattice spanned on two vectors 

il = (Lo,0), /2 = (/21,/22), /22/O. 
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Introducing new variables x = 2-Ky/LQ we reduce the spectral problem for Hb to the form 

(1.2) \ , = E'^\ 



where 



/r, \2(hl\'^ , l^^lmax 

(2.) [-) ,e = h- 



\eB\ . Ifl ■ 
ujc = is the cyclotron frequency , Im = \ is the magnetic length of the system, 

cm V muJc 

iV'Imax = max '^ = TT7\ — y [Lqx / {2t:)) . 

I V I max 

Therefore, the smallness of h means that the characteristic size Lq of the lattice is much greater than the 
magnetic length Im', then e is small if, for example, the electric energy |y|max is comparable with the 
magnetic energy Hujc- 

The smallness of h indicates that the number t] := 022 /h (which is the number of the magnetic flux 
quanta through the elementary cell) is large. 

Such a situation can be realized e. g. in periodic arrays of quantum dots or antidots or in super-lattices [10, 
39,72,93]). 

It is well known that the spectral properties of the operator H depend crucially on the parameter rj. If 
7] = N/M is rational, then the spectrum of H has band structure (in this case, H has the Kadison property, 
see [12, 46]). For each spectral value E of H it is then possible to construct a basis of M generalized 
eigenfunctions {x, q), j = 0, . . . , M — 1, depending on two new parameters q = {qi, ^2), qi G [0, 
Q2 G [0, 1], with the following magneto-Bloch properties (see [12,46,74]): 

(1.3) ^^{x + ai,q) = ^^+\x,q)e-^'''^'^''^'^\ j = 0,...,M -1, 

^^{x + a2,q) = ^^^Hx, g)e-i^{^'i+'^2i/2)^ j = 0, . . . , M - 2, 

Thus the spectral value E becomes a function of {qi, q2) and both E and '^^ depend on the parameters h 
and e (and also on some others); we omit this dependence to simplify the notation. The structure of the 
spectrum of H becomes much more complicate if rj is irrational; in particular. Cantor sets may arise [51]. 

1.2. The Correspondence Principle. The goals and the structure of the paper. We want to exploit the 
small parameter h to obtain asymptotic information about spec H by means of semiclassical approximation 
as /i ^ +0. 

The fact that the parameter 1/ B plays the same role for the operator H as the Planck constant for the 
ordinary Schrodinger equation was first pointed out in [8, 11]. We emphasize that all our assumptions on 
the parameters h, e, and the vectors ai , 02 are essential for our method. If, for example, | Oj | ~ h, then 
V = v{xi/h, X2/h), and instead of "standard" semiclassical methods the Bom-Oppenheimer (adiabatic) 
approximation has to be used in this situation, and it leads to quite different results (see [75]). 

Semiclassical asymptotics are used very widely in problems with discrete spectrum (see e.g. [28,47,61, 
69]) but they are not commonly used in (multidimensional) problems with continuous spectrum. Hence one 
of the goals of this paper is to point out the potential of semiclassical methods for these problems. In partic- 
ular, we want to understand what the conditions (11.31) and (11.41) mean for the semiclassical approximation. 

It is a well known fact that usually there are no universal asymptotic formulas for the spectrum even 
in quite simple situations; one has to describe different parts of the spectrum by different formulas. In 
the discrete case, the various parts of the spectrum arising in this way are referred to as spectral series in 
physics literature; we keep this notation for our situation where the spectrum is continuous. Let us recall that 
the semiclassical approximation realizes the Correspondence Principle: it allows us to describe asymptotic 
properties of the spectrum of the quantum mechanical system via some objects related to the associated 
classical Hamiltonian system. Thus it is natural to ask which "regimes" of the classical phase space should 
correspond to these spectral series; this is the main motivation of this paper. 
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For the magnetic Schrodinger operator (ll.lt . the classical Hamiltonian is 
(1.5) H{p,x,e) = ^{pi + X2)'^ + ^pI + ev{xi,X2). 

For the operator H we want to show that at least in low-dimensional classically integrable situations these 
pre-images are certain subsets or "regimes", to be denoted by Mr, in the phase space which allow a conve- 
nient description in terms of certain graphs. 

Of course, we now have to explain the relationship between the magnetic Schrodinger operator (II. 5t and 
integrable systems, since the Hamiltonian system associated with the latter one is, generally speaking, non- 
integrable. The connection is brought about through the small parameter e: it turns out that the Hamiltonian 
system associated with (II. 5t is "almost integrable", modulo corrections which are exponentially small with 
respect to this parameter. This observation follows from the averaging methods, which were applied in [2, 
16-18, 63,70] to the analysis of the motion of classical particles subject to a strong uniform magnetic and 
certain electric fields from different pints of view. 

The averaging (see|3ll allow us to reduce the original Hamiltonian system, with two degrees of freedom, 
to a system with one degree of freedom i.e. an integrable system. Actually, this reduction does not depend 
on the periodicity of the electric potential v. But in the two-dimensional periodic case averaging leads 
to a Hamiltonian system with phase space the two-dimensional torus T^, and the "reduced" Hamiltonian 
turns out to be a Morse function on T^. This observation leads naturally to a complete classification of 
the classical motion in terms of "regimes" (see sections 4 and 5). The Bohr-Sommerfeld quantization rule 
defines some subsets of M consisting of points and intervals which form the desired spectral series. For 
each point from such a set we can construct a collection of asymptotic eigenfunctions (or quasimodes) of 
the operator H, which are given as power series in the small parameter and are localized in a neighborhood 
of certain domains in the original configuration space (see section 6). 

The next question is to understand how the actual spectrum and the actual (generalized) eigenfunctions 
are related to the constructed spectral series and quasimodes. For instance, in th case of rational flux, the 
constructed quasimodes do not satisfy the Bloch conditions (II. 3t . (I1.4t . and the Bohr-Sommerfeld rule gives 
a discrete subset in contrast to the band structure of spec H in this case. A strict mathematical answer to 
this question is beyond the scope of the "power" approximations used in this paper; we will discuss it only 
heuristically (section 7). In particular, we obtain a heuristic "Weil formula" for the number of subbands 
in each Landau band (section 8) and discuss a connection between our quasimodes and the Harper-type 
difference equations (sections 6 and 8). 

To motivate our considerations, we begin by describing some well known results from the semiclassical 
analysis of one-dimensional periodic Schrodinger operators with a small parameter in front of the second 
derivative (section IJJi. This example allows us already to illustrate the main features of our approach: the 
geometric description of the spectrum by means of Reeb graphs, the semiclassical structure of quasimodes 
and spectral series in problems with continuous spectra, the relationship between different asymptotic for- 
mulas, and the correctness of certain heuristic considerations. 

1.3. Table of notation. We will have to use a somewhat elaborate notation which we summarize here for 
easy reference. 

• e is a small classical parameter in the classical problem, 

• his a. small semiclassical parameter in the quantum problem, 

• is an integer number describing the accuracy of the expansion with respect to e; 

• L is an integer number describing the accuracy of the expansion with respect to h; 

• rj is the number of the magnetic flux quanta through the elementary cell (which we denote by 
r] = N/M in the rational case); 

• ai and a2 are the generators of the lattice F; 

• d = (di, (^2) € is the drift vector of classical trajectories, di/d2 is the rotation number; 

• / = (/i; /2) G is a vector that is conjugate to d, i. e. di/i + ^2/2 = 1; 

• the over-line index~( tilde) indicates a connection with infinite motion; 

• r € N numbers the regimes, (finite motion) and Mr (finite motion), of the classical motion; 

• 9 = (91 1 92) is the vector of quasimomenta; 
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graph 

Figure 1.1. Global classification of the classi- 
cal motion 

• I = {h^h) G is a multi-index indexing closed (contractible) curves on the two-dimensional 
torus belonging to the boundary regimes and implied quasimodes; 

• e Z is the index of open curves or two-dimensional cylinders belonging to the interior regimes 
and implied quasimodes; 

• /i is the quantum number of the Landau level J^^^ ; 

• z/ is the (quantum) number of the "slow drift" action Jg'^^ (it appears in the boundary regimes only); 

• 5 characterizes the neighborhood of the singular manifolds of the classical motions; 

• j E N numbers the magneto-Bloch eigenf unctions; 

• s is the number of the collection of the magneto-Bloch functions, satisfying (11.31) . (11.41) : 

• is the index of a band in the interior regimes. 

1 .4. Averaging, almost integrability, and classification of the classical motion (sections 3-5). The av- 
eraging process gives us an averaged Hamiltonian, such that in new "corrected" symplectic coordinates 
(with generalized momenta Ji, and generalized coordinates <I>, ^2) we can write 

(1.6) H = -K{'ii,^,e) + 0{e~^l'), 

Here Jo{z) is the Bessel function of order zero and Ay = + 9^/9^2- Our main example in this 

paper is connected with the potential 

(1.7) v{x) = Acos xi + B cos{(3x2), 
where A, B, and (3 are positive constants. Then we have 

(1.8) H{3i,]j,e) = Ji + e{AJoiy^i) 005^1 + BJo{(3^/23~i) 005(0^2)) ■ 

Now for almost all 3i, the Hamiltonian H (or "K) may be considered as a Morse function on the two-torus 
= R^/(ai,a2). Using the topological theory of Hamiltonian systems [19,43], for each fixed Ji we 
may separate the motion defined by the averaged Hamiltonian into different topological regimes, which are 
conveniently described by means of its Reeb graph. After a change of the action variable Ji we obtain 
the regimes as the sets of points in phase space which correspond to topologically similar edges of the Reeb 
graph. Then classical motions through points from a fixed regime are topologically similar. It is convenient 
to present the regimes on the half -plane { ( Ji ,E) G ; Ji > 0} where E is the classical energy of the 
averaged system. We give the complete description of the regimes in section 3; the picture for example il.ll 
is given in Fig. II. 11 
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The motion defined by tlie averaged Hamiltonian H takes place in the domain 

(1.9) So = {(ai,^;)GM2;Ji>0, l^-Jil <e(^|Jo(v^)|+5|Jo(/3\/2jT)|)}. 

This domain is the projection of the the actual motion surface, S; any its cutting by the plane Ji = const is 
then homeomorphic to the Reeb graph of the Morse function H (see Fig. II .111 . 
Also, S decomposes into regimes along the curves 

S = Ji ± e|A|Jo(y2j7)| ± 5|Jo(/9 7207)1 1 

which, in this example, form the common boundaries of the boundary and interior regimes. We distinguish 
between the regimes Mr corresponding to finite classical motion and corresponding to infinite classical 
motion. 

Also, it is natural to distinguish between regular and singular boundaries of the regimes, according to 
whether they are external or internal. The internal boundaries may have intersection points which are their 
singularities. 

With each regime, one can associate topological and analytical characteristics. These are the drift vector, 
the Maslov index, the action variables, and the /onn of the Hamiltonian in the action variables. 

In fact, to each inner point of a regime there corresponds a family of closed trajectories on T^, hence a 
family of closed (for boundary regimes) or open trajectories (for inner regimes) on the covering My. To these 
corresponds in turn a family of Lagrangian (or Liouville) tori, for boundary regimes, and Lagrangian (or 
Liouville) cylinders, for interior regimes, in the original phase space To the Lagrangian (or Liouville) 
tori or cylinders (and hence to the regime under consideration) we may associate {a) the vector d = (^1,(^2) 
of the drift in the original configuration space M^, or equivalently the rotation number di/d2 of the related 
closed trajectory on the torus, and {b) the Maslov indices of the related Lagrangian (or Liouville) tori or 
cylinders. 

The rotation number of a boundary regime is equal to 0/0, there is no drift, and there is no preferred 
direction. The Maslov indices of natural cycles on a Liouville torus is equal to 2. 

On the other hand, the rotation number of an inner regime is not trivial, there exists a preferred direction, 
but each cylinder has only one cycle, and hence only one Maslov index, which again is equal to 2. 

Also, in each regime one can introduce a second action variable J2 and find (c) the analytic representation 
of the Hamiltonian in action variables H = ^' (Ji, which depends on the regime. The drift vector 

and the function ^'"(Ji, ^2, changes discontinuously when one passes from one regime to another. The 
correction "K — H does not change neither this rough description of the classical motion nor the general 
asymptotic description of the spectrum, even though a complete analysis of the effected changes may be 
of importance in certain physical problems. In this paper we do not analyze the classical motion in the 
neighborhood of singular boundaries or the behavior of the corresponding part of the spectrum. Thus we 
introduce some small number 6, and remove certain (^-neighborhood of the singular boundary from all 
regimes Mr and Mr. These new sets we also refer to as regimes; we denote them by Mr^s and Mr,5, 
respectively. 

1.5. The global asymptotic structure of the spectrum (section 6). The Bohr-Sommerfeld quantization 
of the regimes Mr^s and M^^^ results in quantized regimes on the "Reeb surface", which after the projection 
onto the energy axis E defines the first approximation in the asymptotics of the spectrum of the original 
operator. The quantization conditions are different for boundary and interior regimes. In both cases, we can 
quantize the variable Ji thus defining the so-called Landau level 

(1.10) 3[^^ = {^ + is)h. 

For boundary regimes, we have in addition a quantization of 32, given by 

(1.11) jH = (i + ^.)/,, 

Here, fi and u are integers with {3[^\3''2^) S Mr^s- However, J2 is not quantized in interior regimes. 
Now consider the numbers J^''^ e), J^*"^) G Mr,s for boundary regimes, and the functions 



6 




Figure 1.3. Global structure of the spectrum 



Hri3\^\J2,e), (ar\a2) G M,. 5, for interior regimes. We then derive the quantized regimes or the spectral 
series on the surface S and their projections onto the domain Sq in the plane {E, 3i)). These sets consist of 
points (for boundary regimes) and intervals (for interior regimes); for example (II. 7t . the result is sketched 
in Fig. O 

Projecting this set onto the £^-axis we obtain the set which describes the first order asymptotics of the 
spectrum of the operator H. Indeed, we have the following result. 

Proposition 1.1. For each r and suitable {/j,, u) or (p,, J2) related to M,.,^ or JAj-^s cmd arbitrary K, L gN, 
there exist numbers 

= H''{3[^\i^\£) + 0{h^ + £^) 

for the boundary regimes and functions 

E^{d2) = H'id[^\32,e) + 0{h'' + e^) 
for the interior regimes, such that the distance between them and the spectrum of the operator H is 0{e^ + 

We have already mentioned that semiclassical methods will also allow us also to construct asymptotic 
(generalized) eigenfunctions (or quasimodes) for the operator H. Actually, each number in the sets just 
described leads to the construction of infinitely many quasimodes, with support localized in a neighborhood 
of the image of the invariant Liouville tori or cylinders in the configuration plane (see Fig. ll.4t . 

Of course, this "degeneration" in the construction of quasimodes stems from the fact that H has contin- 
uous spectrum. We emphasize that the described construction does not depend on the rationality of the flux 
7], i.e. we do not feel any rationality effects. Our results concerning the spectrum of the operator H and its 
quasimodes cannot be improved using semiclassical approximations in powers of the parameters, not taking 
in account tunneling. However, the description of the spectrum on the plane E, Ji by the quantized regimes 
carries more information about the original operator than the description of the spectrum on the energy axis: 
for example, it separates the spectrum according to the different Landau bands, numbered by the index ji, 
and allows us to estimate their width. In the special case (II. 7I) . this width is (see section 6) 

2e(A\U^^)\+B\U(3\f2!^^)\+0{h)). 

1.6. The case of rational flux (sections 7-8). We now turn to the connection between the constructed set 
and the spectrum of H. If e is smaller than h, then the asymptotic Landau bands do not intersect. So in 
this case the constructed semiclassical "asymptotic" spectrum consists of intervals and points on the axis 
E. As in the case with rational flux the spectrum of the operator H has band structure, it means that in this 
situation discrete points define something like "traces" of the (exponentially) small bands, and on the other 
hand there can be (exponentially) small gaps in the intervals in inner regimes, which one cannot catch by 
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Figure 1.4. Projections of trajectories FIGURE 1.5. Separation of subbands 



means of "power" semiclassical asymptotics. Moreover, there exists probably their fuzziness on the surface 
S (and the plane Ji, E) in the direction Ji. 

To clarify this situation (in heuristic level) one can look at these quasimodes from the point of view of 
magneto-Bloch conditions (II. 3t and (ll.4t . It is clear that the described quasimodes do not satisfy these con- 
ditions, but one can use them as a base for constructing the functions satisfying (ll.3t -( ll.4t . The correspond- 
ing pure algebraic procedure (it does not depend on concrete form of the potential v) gives the following 
results. First, it defines certain points on the intervals from the inner regimes describing the "traces" of gaps 
on them. Secondly, it takes off infinite degeneration in such a sence, that for each Bohr-Sommerfeld point 
'Kr{'Ji''\'j'^\£) and quasimomentum q, from the Bloch conditions we obtain M collection of linear inde- 
pendent ("Bloch") quasimodes (It is interesting that the structure of these "Bloch" quasimodes related to 
boundary regimes does not depend on the choice of the coordinates xi, X2. It is not the case for quasimodes 
related to inner regimes: they have the simplest form if the Bloch conditions in coordinates xi, X2 agree 
with the drift vector (rotation number, which is a topological invariant) in such a way, that the latter one is 
(1,0).) On the other hand, the typical degeneration gives the multiplicity M, which means that indeed the 
Bohr-Sommerfeld points corresponds to M exponentially small subbands, separated by exponentially small 
gaps, satisfying to Bloch conditions (ll.3t - (ll.4t . (Recall that M is the denominator of the flux rj.) 

So if one takes a magnifying glass (i.e. construct a more precise approximation) and look at the Reeb 
graph corresponding to a certain fixed Landau level, and its (exponentially) small neighborhood, the follow- 
ing picture appears (see Fig. ll.5t . 

These consideration gives the heuristic "geometrical Weyl" estimates for number 3^ of subbands for the 
fixed Landau level. The idea is that first one has to count them on the edges of the Reeb graph, and then to 
project the result to the energy axis. Final formula for fi-th Landau level in example il.li is J{{3[^^) A^, 
where N is the nominator of flux rj. 

The construction of "Bloch" quasimodes gives also in the first approximation the simple dependence on 
quasimomenta or the dispersion relations. 

At last we obtain difference Harper-like equations, if quantize the averaged Hamiltonian H in naive way. 
It implies the correspondence between constructed "Bloch" quasimodes and quasimodes of the difference 
equations. We discuss this correspondence in sections 6 and 8. 
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Figure 2.1. Typical 
dispersion relations 



Figure 2.2. Asymptotics of the 
dispersion relations 



2. Example: "graph" semiclassical analysis of the spectral Sturm-Liouville problem 

ON THE circle 

2.1. The periodic Sturm-Liouville problem. To explain what kind of results we obtain for the two- 
dimensional Schrodinger operator H, let us consider the spectral problem in L^(M^) with a small parameter 

h>0, 

(2.1) = -fe^ ^ ^1^^ =^^(x), 

where v{x) is a smooth 27r— periodic function. The structure of the spectrum of L is well known, but 
the presence of the semiclassical parameter h introduces certain additional aspects and allows in particular 
to construct certain explicit semiclassical asymptotic formulas for the spectrum. According to the theory 
developed by Floquet, Krein, Gelfand (we refer to the original works [44,58] and to the reviews [42,56,65, 
80,88]) the spectrum of (I2.lt is continuous and, along the energy axis, separated into bands = [E~ , E^] 
and gaps {E^ij^, E'^j^-^), v = 1,2, ■ ■ ■ , E~ < E^ < E~^-^, Eq > mint). (Some gaps may be closed, such 
that £'+ = Each point E S {E~ ,E^) has multiplicity two, and (I2.lt has two linear independent 

Floquet (or Bloch) solutions. It is convenient to parameterize the points in each band by the quasi momentum 
9, < (7 < 1 (under the assumption that one separates the gaps in cases when E^ = E~^^), and to write 
the dispersion relations as 

(2.2) E = E,{q). 

To each point q G [0, 1] corresponds a Bloch function, i. e. a solution q) of (I2.lt with E = E^{q) 

satisfying the Bloch condition 

(2.3) ^^{x + 2tt, q) = e^'^'i^^ix, q). 

Of course, the functions Ey and ^(/(x, q) depend on h, but to simplify the notation we omit this dependence. 
The points q = {) and (7 = 1/2, (g = 1 is identified with q = G) correspond to the ends of the bands and 
give periodic and anti-periodic solutions of (I2.lt . respectively. If E:^ < -EjJ+i then for E = E^ and 
E = E~^^, (I2.lt has only one solution, the Bloch solutions associated to quasimomenta q and 1 — q we. 
complex conjugate. Clearly, '^i,{x, q) ^ L2{^\). Typical dispersion relations are illustrated in Fig. 12.11 

There are no explicit analytic formulas expressing the dispersion relations in terms of the potential v 
even for the simplest potentials (except in the case of finitely many gaps, see [36, Chapter 2], [65, §1.5]). 
Semiclassical asymptotic formulas are available only for large E. We consider now the situation when 
/i ^ +0 in (I2.lt . This problem was studied in many papers and monographs, cf. e. g. [33,35,42,48,56, 
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62,79,91,92], see also [66] and references there-in, but with a somewhat different point of view. We recall 
here several results in a form suitable for us, keeping in mind their multidimensional generalization. 

2.2. The asymptotics of the spectrum. To simplify the discussion we now assume that v is analytic and has 
only one non-degenerate minimum point Xmin on S^; we may and will assume Vmin '■= ^^(^^min) = 0. Then 
there also exists only one global maximum point Xmax of v; we put Wmax = ^^(a^max)- Under this assumption 
the dispersion picture is divided in four domains. The bands situated under timax become exponentially 
small with respect to h, and the corresponding dispersion curves are almost horizontal segments Ey = const 
mod 0{h°°) (see Fig. l2.2l) with distance 0{h) between them; this means that the number of bands increases 
as h tends to zero, and u is allowed to be large. Let 6 denote a small number independent of h. 

Proposition 2.1. (a) For Ey{q) < v,„ax — S we have 

(2.4) E,iq)=Ei^,+oih), 
where Ei^y is defined by the Bohr-Sommerfeld rule 

1 r+ I 1 

(2.5) - / \ Ei^y-v{x)dx = h{- + u), 



7rJ,_ V - ^ ^ ^2 

with x±(Ei^i^) solutions of the equation v{x) = Ei^^ (see Fis. 12. 3 1 . 

If uh is small enough then Ei^y is also small and one can simplify (12.51) using a Taylor expansion; this 
leads to the "harmonic " oscillator approximation for Ei y, 



(2.6) ^1,^ = +l/)cJo + 0(/l^), wo = \/lv"(Xn,in)- 

(b) Let v,nin + 5 < Ei^y < v^ax - ^, then 

(2.7) Ey{q) - E- = cos(27rg) + l)e-''/^(l + 0{h)). 



Here uj{E) = 2Tr{f^'^ --j===dx) ^ is the frequency {see subsection 12. 4t . and p = j'^^ ^ v{x) — Ei^i^dx 

is the Agmon distance [1]. 
Ifvh is small enough then 

24-+5/2^0-+3/2/jl/2-. 



\/v(x) 



sm I 



{{-ir+^ cos(27rg) + 1)6-"/'^ (1 + 0{h^/^)) , 

with Agmon distance p = ^"""+2'^ \/v(x) — Ei ydx. 

For the proof we refer to [91], see also review of results in [66]. 
Remarks. 

(1) If ~ 1/h then o{h) in (I2.4t can be replaced by 0{h?); the estimate o{h) appears during the passage 
from "small" to large "large" u, see [56,66]. 

(2) In these asymptotic formulas the potential v appears only through the frequency and the Agmon 
distance: the dependence on the quasimomentum q is the same for different potentials. 

(3) If we formally take the limit uh Oin (I2.7t for uh « 1, we obtain (I2.8t . but these arguments are 
not rigorous, so (I2.8t has to be proved by other means. In such a situation we say that formula (I2.7t allows 
a formal limit. 

(4) The subtraction of 1/ sin ( ) from the integrand in (12.81) is just one possible type of regularization. 

In the upper domain, on the other hand, the gaps above i^max become exponentially small, the bands have 
length 0{h) and almost cover the spectral half axis. 
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Proposition 2.2. For Ey{q) > Vmax + S we have E'^ — E^j^^ = Oy{h°°), and the ends of the gaps are 
defined by a Bohr-Sommerfeld quantization rule in the form 



1 /■^'^ hu 



(2.9) E+ = 8,, + 0,ih^), — ^l, - v{x)dx . 
We also have the following asymptotic formulas for the dispersion relation in the v-th band (Fig J2.2t : 

(2.10) E,{q)=E2,u{q)+0,{h^), ^ ^ E2M - v{x)dx = I^{q,h), 
where for even v = 2s 



(2.11a) h{q.h)-- 
and for odd u = 2s + I 
(2.11b) I,iq,h) 



h{^+q), 0<g<i, 



h{i^ + q), \<q<l. 

For the proof of ( I2.9I ). see [35, Appendix B], or [42, §7]. Formulas ( 12. 11 at and (12.1 Ibt in a slightly 
different form are contained in [30] or [56, Part 11]. 



Remark. The description of the splitting E~_^-^ — E^ between the ends of the gaps is not so simple in this case 
as in (I2.7t . In the physics literature, this splitting is associated with the so-called "over barrier" reflection. 
Under additional assumptions [35,42,91,92], the splitting has the more explicit form: 

(2.12) E;^, - Et = ^^&!^e-Pl^{l + 0{h)). 

vr 

The definition of the Agmon distance p, however, is now quite different. 

Of course, there is a transient layer in a neighborhood of Umax where the band length is comparable with 
the gap length; formulas (12. 5t . (12.71) . (I2.10t . and (12. lit are not valid there. We do not consider this situation 
(see, nevertheless, [15,66,92]). 

2.3. Quasimodes and Bloch solutions. The behavior of the Bloch solutions differs sharply in the lower 
and upper domains. Moreover, it is natural to isolate a certain neighborhood of the bottom of the lower 
domain, because the behavior of the corresponding eigenfunctions is also different there. The asymptotic 
formulas depend, of course, on the accuracy of the approximation: they must be more complicated e. g. in 
case of the subtle dispersion relations (I2.8t . M.9\ . 

Recall the following definitions (see e. g. [49-51, 61, 67-69]). 

Definition 2.1. Let L > 1 be a real number. 

• A pair E^) is called sl formal asymptotic solution or quasimode of order L, relative to some 
function space J, if 

(2.13) ||(L-^^)^'^||:r = 0(/i^). 
We can use, for example, J = C(M), or 3~ = L'^{M.x)- 

• Let be a solution of the equation (L — E')^ = 0. The function ^'^ is called an asymptotic part of 
order L of ^' if ||* - ^-^Hgr = 0{h^) as /i ^ 0. 

• Let 1 1 ^^11 > c > as /i — > 0. The function is called a leading term of the quasimode if 
||^i_^0||^^o(l)as/i^O. 

Remarks. 

(1) Note that in the definition of quasimode, the Bloch condition (12.31) is not required. 

(2) Definition (I2.lt describes so-called "power" or "additive" asymptotics; these notions are used in 
contrast to "multiplicative" asymptotics, which we will define later. 

(3) An asymptotic part contains more information about the true solution of (I2.lt than a quasimode, even 
though both concepts can coincide in specific examples. Nevertheless, one can derive information about the 
spectrum of L from quasimodes. The following proposition is essentially well known. 
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Proposition 2.3. Let be a smooth function and G M with the property 

is a quasimode of L of order L in L^(M), and > c > as h ^ 0; 

or 

(b) (^'■^, E^) is a quasimode of L of order L in vr, vr], ^'^ satisfies (I2.3t . and ||^''^||L2[-7r,7r] > c > 
Oash^O. 

Then the distance between E^ and the spectrum of L is 0{h^). 

Proof. The proof is well known for the case (a), see, for example [67, Lemma 1.3] or [69, Lemma 13.1]. 
Let us give the proof for the case (b), which is a simple generalization and probably also known. 

For any function Lp satisfying (I2.3t we have Hv'llLaj.nj^^^] = \/^ll'/'llL2[_7r,7r]> so this holds, in par- 
ticular, for for {^^y , and for the discrepancy f := {L — E^)'^^. Let us choose a smooth cut off 
function e with < e < 1, e(x) = 1 for x G (— 7r,7r), and e{x) = for x ^ (— 27r,27r), and let 
|e(x)| + |e'(x)| + |e"(x)| < ci for x G M. For m G N we put em{x) := e{x/m). Since L is self-adjoint, 
we have 

(2.14) ||e^^'^|U2(K)dist(specL,£;^) < ||(L - ^^)(e„^'^)||i2(K), 

and the left-hand side satisfies 

(2-15) ||em^'^||i2(R) > ci/m, 

by assumption. To estimate the right-hand side we use 

Then 
(2.16) 

\\{L - ii;^)(e^M/^)||i2(K) < ^/2^||/||i2[_,,,] + 



m\/m ' ' Wm 



dx 



Combining (12. 14t . (I2.15t . and (12. 16t we obtain in the limit m — > +oo the desired estimate 

dist(spec L,i?^) < ^||(L-ii;^)M/^||i2[_,,,] =0(/i^). 

□ 

To describe the asymptotics of Bloch functions let us start from the simplest level of complexity related 
to ( 12.61 ). One obtains the following picture: Bloch functions associated to the lowest bands are localized in 
0(\//i) -neighborhoods of the minimum points x^in + 27r/, / G Z, of the potential v, where they coincide to 
first order with the eigenfunctions of a harmonic oscillator. More precisely, in some 0(\//i) -neighborhood 
of Xmin one has the following formula for the leading term in the asymptotics of all Bloch solutions: 

(2.17) (3;) = exp ( — )}ly{ 



where C is a normalizing constant and H,^ denotes the u-th Hermite polynomial, whereas the Bloch func- 
tions are 0{h°°) in all other points of the segment [xmaxi ^^max + 2vr]. This together with the Bloch condition 
completely defines a leading term in suitable neighborhoods of all other minimum points Xmin + 2ttI, I G Z, 
by the formula 

(2.18) ^o'(x,q) = J^e2-«Vo"(^ - 2^0- 

leF 

More precisely, for any Bloch function ^''^ there exists a function ip'^ such that 

00 

(2.19) ^'''(x,g) = ^ e2'^*«V(a;-2vr/); 

/=— 00 

this is the so-called Gelfand representation (see [44, 88], [80, XIII. 16]. Let us record the fact that tp'^ = 



12 




Figure 2.3. Structure of Bloch functions 



Using the terminology introduced above, one can prove that ( 12.171 gives asymptotics of certain quasi- 
modes of order L for (I2.lt . and the functions (I2.18t are the leading terms of asymptotics of the Bloch 
solutions. In a way, (12. 18t presents the asymptotics of modes via quasimodes, and the approximation (12. 17t 
and (12. 18t allow us to derive (I2.6t . Note that (12. IBt gives more information about the Bloch solutions 
than (12. 18t but no better spectral information than (12.61) . 

The Bloch solutions corresponding to the higher bands are localized in a neighborhood of the segments 
[x_ + 27r/,x+ + 2-kI], I G Z, where x± are solutions of the equation v{x) = Ei^y introduced above; they 
can be represented in the form ( 12. ISt . This means precisely that in inner points of the interval (x., 
a leading term of all Bloch solutions is given by 

(2.20) ro{x):=- ^^-T7l(cos(- / J E,,, - v{x)dx + -) + 0{h)) , 

with C a normalizing constant. 

In a neighborhood of the turning points x_ and X4., the functions 4'q{x, h) are given in terms of Airy 
functions and have large amplitudes; but they are still 0{h°°) outside certain neighborhoods of the segment 
Thus it follows from ( I2.18t again that there exist gaps in the asymptotic support of the Bloch 
solutions (see below for the definition). A global uniform "power" asymptotic of ip^ can be given in terms 
of Maslov's canonical operator (we will return to this representation later). Using quasimodes as before, 
one derives the spectral information given in ilAl and (I2.5I) . 

It is convenient to use some terminology taken from the theory of short-wavelength approximation 
in optics Let us consider a certain asymptotic solution ^'(x, /i). The closure of the domain where 
lim/i^o "^{x, /i) 7^ is called its asymptotic support or light region. The domain where "^{x) = 0{h°°) 
as /i — > is called the shadow region. In some neighborhood of the boundary (this neighborhood is small 
together with h) ^'(x, h) has order h^; sometimes this neighborhood is called the penumbra (this defini- 
tion, of course, is not rigorous). So for (12. 17t (or (I2.18t ) the light region is the union of the minimum 
points Xmin + 27r/, / G Z, all other points belong to the shadow region, and the penumbra is some neigh- 
borhood of {xmin + 27r/, / G Z}. The light region for the asymptotic solutions related to the higher bands 
consists of the union of the segments [x^ + 2ttI, x+ + 27r/], / G Z, all other points belong to the shadow 
region, and the penumbra is the union of certain neighborhoods of the turning points x_ + 27r/, x+ + 2^1, 
Z G Z. In quantum mechanics, the shadow region is also sometimes called the under-barrier region. 

Now let discuss the representation (12.181) in greater detail. 

Proposition 2.4. We fix e > and denote by e some smooth cut off function with e{x) = 1 for x G 
{b — e,b + 2tt + e) and e{x) = Ofor x ^ {b — 2e,b + 2tt -\- 2e) {the number b is defined later). 

(a) Let v £ N be a fixed number, then the Bloch function associated with the u-th band has the form 
( I2.19L where ip'^{x) coincides up to 0{^/li) with the function (12.171) in a certain neighborhood of Xmin- 
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Outside this neighborhood, in the interval (b — 2e, 6 + 27r + 2e) we have 



(2.21) 



/ 95!^+l,,,i^+l 1 /tnn f ^"^""" ! 



Here ^{x) := \ . ^/v{x)dx\, ^(x) = ^'{x) = ^/v(x) if x > x^in and ^{x) = ^'{x) = — y^f (x) 
otherwise. The point b is defined as the unique solution of the equation ^{b) = ^{b + 27r), /. e. 
v^dx = £+f ^)dx. 
(b) Let ci < hv < C2 for some ci, C2 > 0, G N. Then the Bloch function associated with 
the vth band has the form (I2.19I ). where ip'^{x) coincides up to o{h) with (I2.20t in the interior of the interval 
(x_ , Moreover, 



(2.22) ^P^ix) 



V^{x)-E,,,dx) (1 + 0(/.)) , x<x.-e, 



and b is defined by the equation v{x) — Ei^i, dx = J^^ ^ v{x) — E\^y dx. 

The normalizing constants C in (I2.21t and (12.221) are the same as in (12.171) and (12.201) . 

Sof/i asymptotics (12.211) a;«(i (12.221) admit differentiation with respect to x, i.e. in these formulas 
{0{h))l = 0{h). 

For the proof, we refer to [42, §§ IV. 1, IV.4]. 
Remarks. 



(1) Formula (12.211) does not admit a formal limit for x Xmin', in particular, (12.211) contains only the 
highest degree term (y^ljo/2(x — Xmin)) of the corresponding Hermite polynomial, but the other terms 
also play a role in a neighborhood of the minimal points. Of course, it is possible to extend the construction 
accordingly, but, as we will see below, it is not necessary for obtaining the dispersion relation il.Sl . The 
presence of the tan-like term in (12.2 It is only one possible way of regularization; another way of regular- 
ization has been used in [48]. 

(2) One has different constructions for the Bloch functions corresponding to the bottom lower and the in- 
ner lower domains. In the first case, they are defined by a real-valued phase and decay exponentially with 
h, while in the second case one needs complex phases, i. e., in this case the Bloch functions have both 
oscillating and exponentially decaying parts. This phenomenon reflects deep properties of the asymptotics 
given by Maslov's canonical operator; this distinction appears more clearly in multidimensional problems. 

It is necessary to emphasize that it is not complicated to obtain the asymptotic formulas for the spec- 
trum in this one-dimensional situation, but it is difficult to prove the formulas for the true asymptotics of 
the Bloch functions. The standard way of doing this in the one-dimensional case is based on WKB meth- 
ods for ordinary differential equations and matching solutions in the complex plane, see [42,79,89,91,92]. 
These methods are applicable in both bottom and inner parts of the lower domain, and allow to obtain also 
the corresponding dispersion relations. But up to now there is no rigorous generalization of this method to 
multidimensional problems. On the other hand, there are some methods (see e.g. [1,33,48-51,67,68,84,85]) 
which are applicable also to multidimensional spectral problems (like tunneling problems or problems with 
purely imaginary phase), but they work in the bottom parts of the spectrum only. One may call all these 
methods "semiclassical approximations" (although not in sense of [60]), because they use certain objects 
from classical mechanics. 

(3) Usually, Bloch functions corresponding to the same band and to the quasimomenta q and 1 — g are 
normalized in such a way that their Wronskian is equal to 2i. This leads to normalizing constants in (I2.17t 
and (12.201) which are exponentially large in h. Otherwise, the behavior of the Bloch functions is quite 
strange: in each segment [2ttI, 27r(/ + 1)] some their linear combination is 0{hP°). The appearance of large 
normalizing constants destroys this effect. 
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Clearly, the exponential smallness of the lower bands makes it difficult to calculate the spectrum and 
the Bloch functions numerically. 

(4) If for some solution of (12.11) we have a representation (l + 0(/i))^'o which can be differ- 
entiated in X, then the function ^I'o is sometimes called a multiplicative asymptotic of ^; both formulas 
(I2.21t . (I2.22t together with (I2.19t provide examples. In contrast to additive asymptotics, multiplicative 
asymptotics make sense also in the shadow region. Multiplicative asymptotics are sometimes also called ex- 
ponential or tunneling asymptotics, because knowing them allows to construct asymptotics of the spectrum 
with an error 0{e~'-^/^) which is necessary to deal with tunneling effects. 



Let us show now that from the formulas (12.211) and ( 12.221) it is easy to derive the dispersion rela- 
tions (123 and (|23. Of course, this can be done using matching solutions in the complex plane and the 
one-dimensional WKB method - as mentioned above, - but we give here a derivation based on the simple 
integral formula suggested by I. M. Lifshits (see [62, §VI.55, Problem 3]). 

Recall that if {^i,Ei) and (^'2,-^2) are solutions of dTTl and ^'i^2da; / 0, then 



(2.23) 



i?i — Eo 



Ev{qi), i = 1,2, a 



b, P = b + 2Tr. Note that both 



Let us choose in (I2.23t 'l'i(x) = ^"{x, qi), Ei -- 
^1 and ^2 are defined by ( 12. 19t . and thus we have 

(*i^'2 - ^2^'i) = 2((^'^)'(6 + 27r)V'''(6) - iJ^ib + 2^)(V'^)'(6)) (cos(2^gi) - cos(2^gi)) , 

because the supports of the functions ip'^{- — 2-kI) do not intersect. In the inner lower domain, for the de- 
nominator of (I2.23t we have 

h27r _^ _ />fe+27r 

^'i(x)^'2(x) dx = 



(2.24) 



and, therefore. 



h,;2=-i,o,i 



2TTi{liqi+l2q2) 



i^^ix - 2-nli)'4)'' {x - 27r/2) dx 



['"^^\r{x)fdx + 0{h^) = t^^\rQ{x)fdx + o{h), 

Jb Jb 



(2.25) 



E^iqi) - EM = 2/i' 



{^j'^Yib + 27T)ij''{b) - V{h + 2^){'4,^)' {b) 



t'WoYdx 
(cos(27r(7i) — cos(27r(72)) (l + . 



For the ground lower domain we have ( 12.241 and ( 12.251 ) with o{h) replaced by 0{Vh). 
To prove (12.71) and (12. 8t one has only to substitute (I2.21t . M.221 into (12.251) . 

We observe again that the main term of the dispersion relation is independent of v; the potential appears 



only in its coefficients in terms of higher order in h. To determine the nominator in (I2.25t . one should use 
a multiplicative asymptotics, while the denominator is defined by power asymptotics. 

The Bloch functions for all quasimomenta in the bands from the upper domain are bounded as /i ^ and 
oscillate everywhere on M^, there are no gaps in their asymptotic support; they can be expressed by means 
of simple formulas outside neighborhoods of size 0{h°°) of the points q 



0,i: 



(2.26) 



la V^-^Al^h) -v{x)dx) 



+ 0(/i) 



, {E2Aq,h)-vix)f' y 

HereC"''=^(g) and a are normalizing constants, -£^2,1/ (9) is defined by ( 12.101 ) and ( 12.111 ). and one has to take 
signs -I- and - according to g G (i, 1) and g G (0, ^), respectively. 

The formula (12.261) does not give asymptotics of the solutions of (I2.lt in the points q 
the ends of the bands. Due to resonances and tunnel effects between these points, the asymptotics of the 
true eigenfunctions (the periodic and antiperiodic solutions) are given by the even and odd combinations 



of the functions (I2.26t provided that the constants and a in (I2.26t are chosen appropriately, see [35] 



0, ^, i. e. in 
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for details. In some cases, the constants C' can be expressed through each other; to do this one can 
normalize the corresponding Bloch function "^'^{x, q) by the condition ^''^(0, q) = I (see [36]). 

So for these eigenvalues, (12.261) defines quasimodes but not asymptotics of the modes. 

The Bloch functions related to the transient domain are close to the latter ones, but in a neighborhood of 
the critical points a^^ax + 27r/, I G Z, one can express them by means of certain special functions; e. g. if 
a^max is a non-degenerate critical point of v, then these are the Weber (parabolic cylinder) functions. 

The results we have mentioned so far are obtained by a variety of techniques but with different lev- 
els of complexity. Thus it is considerably more difficult to derive the dispersion relations (12.71) and (12. 8t 
with exponentially small bands — using "multiplicative asymptotics" corresponding to tunnel effects — 
than (12.41) and (12. 6t . The analysis of the transient domain — which we have not explained here — becomes 
even more complicated. 

Some of the methods mentioned above have been extended to problems in higher dimensions but not 
in a systematic way. For such an approach, from the general philosophy of quantum mechanics we should 
expect a correspondence between certain characteristic parts of the spectrum of L (so-called spectral se- 
ries) and certain characteristic geometric objects in the phase space of the classical motion. In our one- 
dimensional example, the classical motion is integrable, such that inspiration gained here can be expected 
to extend at least to the generic integrable case, and that is what we want to explain. 

In the case at hand, the spectrum of L may be decomposed into four domains having similar asymptotic 
behavior as detailed above; these are spectral series. We are now going to show that the presence of different 
types of asymptotics naturally corresponds to a decomposition of the phase space into "regimes" which each 
allow a simultaneous treatment of the flow. This decomposition, in turn, is characterized by a single graph 
which, in this example, coincides with the Reeb graph of the corresponding classical Hamiltonian. 

2.4. The graph of the classical motion. We now want to construct classical preimages of the spectral 
series described above. To do so, we give a suitable classification of the classical motion and establish 
a relationship with the "quantum motion" defined by (I2.lt . Thus we have to consider the corresponding 
classical problem defined by the one-dimensional Hamiltonian 

(2.27) 

The related Hamiltonian system 

P 

(2.28) 

X 

can be considered from two points of view: (1) as a system with phase space M.p^^; (2) as a system with 
phase space the cylinder ^ :=Rp x Sl, such that 3, is the universal covering of Qp^^. 
Then we can distinguish the following types of the trajectories: 

a) closed trajectories on which correspond to closed contractible trajectories on Qp ^ (on these 

we have Vmin < H < Vn^^x); 

b) open trajectories on M^^ which correspond to closed but not contractible trajectories on Qp ,^ (on 
these we have H > fmax); 

c) the stable minimum points (0, x^i^ + 2ttI) on or o^i Qp^x'^ 

d) the saddle points (0, Xmax + 27r/) on M^^ or on Qp ^ and the singular manifolds (separatrices) on 
Mp J. or Qp x, which belong to the "singular" energy level fmax- 

This correspondence can be easily illustrated if one imagines that the trajectories of the Hamiltonian 
system are level curves of the height function of the deformed cylinder, see Fig. 12.41 

Both phase pictures decompose qualitatively into the stationary point(s), the separatrix, and the three con- 
nected components of the complement of their union; obviously, two of the three components are equivalent 
under the map p 1— > —p. Relating this to the energy function, we see that the stationary point(s) corre- 
spond(s) to Viain while the separatrix corresponds to Wmax- Hence the Reeb graph G [19] of H describes 
the situation nicely, cf. Fig. 12.41 The Reeb graph is constructed as follows: each connected component of 
the level set of H corresponds to a point of this graph, and connectivity in this set is introduced in a natural 
way. In our case, the Reeb graph has four vertices, Vmin, Vmax, and 002/3, and edges ii := (wmin, I'max), 



H{p, x) := + v{x). 



-v'{x) 



2p 



x = —2v'{x), 




Figure 2.4. Reeb graph and trajectories of Hamiltonian system 

^2/3 •= (^'max, 002/3)- Morcover, each edge may be identified with an interval of the energy axis. We ob- 
serve next that for each edge of the graph we have to introduce a different action variable which we denote 
by F, where r = 1, 2, 3 numbers the corresponding edge. For the edge ii we have 

(2.29) l\H) = ^ I pdx = - I ^ - V{x)dx, H < Vraax- 

27r JaI tt J^_ 

For ^2/3 we have: 

(2.30) 1^/3 (if) = i- J^^ - V{x)dx, H > Vrr^. 

Introduce the action variable in the saddle points: 

/^"^(^^max) = lim I\H) = - \/ Vmax - v{x)dx, 

-ff-»Umax-0 TT Jq 

(^^max) = lim = — / V^max - v{x)dx. 

Obviously, limjy^„„,,+o /H^^) = 0, and \iuiH^+ooI^'^{H) = +cx), so G [0, /i+(t;max)] and /2/3 g 
[j2+/3+^ +00). One has the "Kirchoff law" /^+(?;max) = ^^"(^^max) + ^^"(t'max), such that /^+(^;max) = 

2/2-/3- 

Since the functions F are continuous and strictly increasing, we can invert them to find the dependence 
of H on I for each edge, 

(2.31) H = ?{^(/), 

where 'K'^{I) = ^^{I). Next we will also parameterize the trajectories by the action variables, separately 
for each edge i^, r = 1, 2, 3. 

We have obtained three open subsets in the phase space corresponding to the edges v to be denoted 
by Mr, r = 1,2, 3; these are the regimes mentioned above. 

In Ml, we have only closed trajectories grouped together by their images in Q^. These curves can be 
parameterized by the action variable as follows: 

Ajil) = {pu{l,t),xuil,t)), I e (0,/!+), t G [0,T(/)], I e Z, 
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where [pii{I, t),xii{I, t)) is the solution of (I2.28I) satisfying 

pi/(I,0)=0, xio(/,0) =x_ (Jf !(/)), 

and 

xi/(/,0) =xio(I,0) +27r/. 
Thus all the trajectories A^^ are uniquely determined and periodic with period 

P+(:k1(/)) . 
Til) = / , dx; 

all K] cover the same trajectory on Q^. 

In M2 (and hkewise in M3) we obtain quite similarly families of trajectories 

= {p2{I,t),X2(,I,t)), IG (/2-,oo), f GM, 

where {p2(I, t),X2{I, t)) is the solution of (I2.28t satisfying 

P2/3iI, 0) = ±^?{2/3(/)-t,^,„ X2/3(/, 0) = Xn,ax. 

A2/3 covers a unique trajectory, A^/^, on and enjoys the periodicity property 

P2/3 {I, t + Til)) = P2/3 {I, t) , X2/3 (/, t + r(/)) = X2/3 (/, t) ± 27r, 

where now 

1 r 1 

T(/) = - / = dx. 

All the trajectories A^^, A^/'^ are one-dimensional Lagrangian manifolds. The closed curves A^^ have 
a Maslov index equal to 2; the curves A^/^ are open and their Maslov index is not defined. 

Finally, we see that the phase space is separated into regimes corresponding to edges of the Reeb graph 
for H. Trajectories from the same regime have similar topological characteristics. Singular manifolds form 
the boundaries of the regimes. 

Of course, for a general potential v (even required to be a Morse function) the Reeb graph can become 
very complicated. It is impossible to give a "generic" description of the Reeb graph because there exists no 
"generic" potential. But obviously the procedure described above is applicable in any case. 

2.5. The relationship between the graph and spectral asymptotics. Now it is quite easy to see that our 
regimes are suitable objects for semiclassical quantization or, more precisely, that they explain the spectral 
series of the Sturm-Liouville problem (I2.lt as described in 12.21 above, corresponding to the four energy 
domains. Indeed, we set up the following relationship (Fig. l2.6t : 

bottom lower domain < — > bottom part of the regime Mi ; 

inner lower domain < — > inner part of the regime Mi ; 

"transient" layer < — > some small neighborhood of the boundary between Mi, 

and M2/3; 

upper domain < — > the regimes M2 and M3. 

Keeping in mind the previous explanation concerning the "transient" layer let us introduce new regimes 
Mi_5, Mr,5, r = 2, 3 which are the "old" regimes but without certain (5-neighborhoods of the singular points 
/ = I^"*", / = /^~/^~; we will describe the semiclassical quantization in these domains. 

Consider first the regime Mi^^ (related to the edge ii). The Bohr-Sommerfeld rule (I2.5t in this situation 
may be rewritten in the form 

(2.32) I = &~^ = h{^ + u^ 

and gives the "quantized regime" or the spectral series corresponding the the regime Mi 5 (or to the edge 
ii). The non-negative integers v are chosen in such a way that /('^^ G [0,/^^ — 5\. Hence v ~ if 
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Figure 2.6. Relationship between the Reeb graph 
and spectral series 



I > 5 > 0. The map (|233 together with the closed curves A,^ gives the set Ei^^ = IK'" of 
"asymptotic" eigenvalues (12.51) and the quasimodes for the operator L; in fact, one can prove that for each 
L, Z G N one can find the numbers (independent of /) 

(2.33) ^ = ^1,^ + 0{h^) 



and a family of quasimodes ^^'^{x, h) of (I2.lt of order L such that 



(2.34) 



where tTx denotes the projection onto the x-plane^ This general construction is well known (see e.g. [42, 
56,69]) and may be carried out using Maslov's canonical operator on the curve A^i (for 

one can use the real canonical operator, in case I*^'^) ^ it is necessary to use the complex 
canonical operator). For the leading term one has 



(2.35) 



i:'{{x,h) 



If I^'^^ is small, then the last formula for / = is (I2.17t . If I^'^^ > k > 0, then outside some small 
neighborhood of the intervals [x_ + 27r/,x+ + 2ttI] = TTxAj{I^'^^) we have ^/;]^(x,/i) = 0{h°^), and for 
ipQix, h) one has formula (I2.19t in the inner points of the interval (x__ , x+). 
One may also express ip\{x, h) via V'o (a^j h) by the formula 

(2.36) 'ilj];{x,h) = i;^{x -27rl,h). 

Now let us return to quasimodes and the Bloch conditions (I2.3t . We know that (I2.18t holds in the case 
at hand, but its proof uses additional nontrivial asymptotic considerations, and some of them are not yet 
available in multidimensional situations. But let us give some simple heuristic and purely algebraic argument 
how to obtain (I2.18t with the ansatz 

oo 

(2.37) %ix,q,h) = Ciiq,h)roix-27rl,h), 

l=—oo 

where Ci{q, h) are unknown coefficients. Requiring the Bloch condition we find 

oo oo 

Ci{q)ro{x-27r{l-l),h)=e^^"^ ^ QiqWoix - 27rl,h). 

/=— oo l=~oo 

27rl,h))^ 



If the system (^TpQ {x 
(2.38) 



has suitable basis properties, we conclude 

Ci+^iq) = Ci{q)e^^"^, 



^Eq. <2.34t means that ■i/)p^(a;) ^ as ft — > for any x ^ tv^A] . 
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Figure 2.7. Exact dispersion relations 

hence C; = e^'^*'?'C^, where is a normalizing constant and q G [0, 1), and we obtain formulas ( 12.181 ). 
(12. 19t . and as corollary the structure of the dispersion relation (12.251) . This consideration does not depend 
on L, the degree of approximation. 

So in this case the used semiclassical method gives a 0{h°°) -approximation of the dispersion relations 
and the asymptotics for the Bloch functions. 

Now consider the regimes Jiir^s corresponding to the edges v, r = 2, 3. There are no cycles on A^'^(I), 
and for each / G M,. 5 and arbitrary large L one can write the following formula for the asymptotic solutions: 



(2.39) 



^p^{x,h,E) = & 



±1 la - V{x)dx 



+ 0{h) 



where the sign + corresponds to M2, the sign - to M3, and a and are some constants. The function if}^ 
is associated with the spectral value 

(2.40) E^{I, h) = ■K^{I) + 0{h^) = 'K^{I) + 0{h^). 

Requiring now the Bloch condition for the functions ip^, we derive the dependence of / on g as 

(2.41) l'^{q,h) = h{n + q), I^{q,h) = h{n - q), n G Z. 
This dependence also implies the dependence of the energy on the quasimomenta, 

(2.42) E^/^'^ (q) = E^ (/^/^ {q,h),h). 

Recall that the points of the spectrum corresponding to periodic and anti-periodic solutions of (I2.lt lie on 
the ends of the bands. Applying this fact to the function ( I2.39t one immediately obtains the points 

(2.43) I^"^ = hv/2, uel 

from Mr,5, r = 2, 3; the corresponding energy levels E^(^I^'^^) are therefore 0(/i°°)-approximations of the 
gaps. 



Combining now (12.4 It . (I2.42t . and (I2.43t . we arrive at the dispersion relations (I2.10t - (I2.12I) . 



Note that points (I2.43t with even u may be obtained by means of the Bohr-Sommerfeld quantization 
of the non-contractible closed preimages of A^'^(/) on the cylinder Qp.a- This fact has a rather simple 
explanation: these points correspond to periodic solutions of (I2.lt . and only these Bloch solutions descend 
to functions on the cylinder Qp ^- Anti -periodic solutions do not descend to functions on Qp ^, but only to 
the enlarged cylinder Qp ^ = Rp x S^, x ^ [0, 47r). The Bohr-Sommerfeld quantization rule on Qp^ then 
gives exactly the points (12.431) . 

Figure 12.61 shows the relationship between action variables, quasimomenta and energy. This picture 
together with formulas (I2.4t - (l2.6t . (I2.9t - (I2.1 It . (I2.18t . (I2.25t contains the maximal information about the 
spectrum of L which can be derived from additive asymptotics. 

The precise structure of the dispersion relations is sketched in Fig. 12.71 it is not accessible in details by 
these methods. 
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2.6. The Weil formula. To conclude this section, let us suggest a heuristic method for calculating the 
number N{E) of bands on the half-line (—00, E) (the Weyl formula). 

\f E < Umin, there are no real trajectories of the Hamiltonian H{p,x), no points on the graph G and 
N{E) = 0. 

If G (I'mini^'max " S), then the number of bands approximately coincides with the number of the 
Bohr-Sommerfeld points I^, (12.321) in the interval (wmim E), i.e. with I{E)/h. 

\i E > fmax + 5, then we have gaps on the edges 12 and ^3, but by symmetry their projections to 
the energy axis coincide and one has to take into account only one edge, say 12, which gives N{E) = 
I^+/h + 2{I{E) - fl-)/h = 2I{E)/h. 

Last two formulas have common geometric interpretation: hN{E) is approximately equal to the square 
of the set < X < 27r, H{p, x) < E, i.e. the set covered by the trajectories of (I2.28t with energy not greater 
than E. 



3. Classical averaging 

Now we return to the spectral problem of the magnetic Schrodinger operator (ll.lt . We will use ideas 
closed to those collected in the previous section, but we will start directly with the classical problem. First 
we want to show that the presence of the small parameter e renders the almost integrable system. Basing on 
this fact, we give a global geometric classification of the classical motion in the following sections. 

Consider the classical problem in the phase space R^^ induced by the operator H and defined by the 
Hamiltonian (11.51) : 

(3.1) H = Ho + ev{xi,X2), Hq = ^{pi + X2f + ^pI 

The projections of the trajectories of the Hamiltonian system with free Hamiltonian Hq onto the (xi, X2)- 
plane are the cyclotron circles, see e.g. [2, 16, 18,63,78], and they induce new canonical variables in the 
phase space: generalized momenta /i, yi (or P, yi) and generalized positions Lpi, y2 or (Q, 2/2): 

xi = Q + yi, pi = -y2, x2 = P + y2, P2 = -Q, 

(3.2) , — , — 

P = y 2I1 cos , Q = V 2/1 sin Lp\ , 

such that 

dp\ A dx\ + dp2 A dx2 = dl\ A dc^i + dy\ A dy2 = dP A dQ + dy\ A dy2- 
The variables P, Q (or /i, 99) describe fast rotating motion around slow guiding center with coordinates yi, 

y2 [63]. 

In these variables, the Hamiltonian H takes the form 

(3.3) H = Ii + ev{y/2hsmipi + yi, \/27i"cos(^i + 7/2), 

and furnishes probably the simplest example where the averaging methods (see e.g. [4, 16, 17,71]) can be 
successfully applied. The averaging procedure for the Hamiltonian H was first applied by van Alfven [2] ; 
later, it was used in numerous works (usually, not in the variables (/i, ipi,y), y = (yi, y2)> see e.g. [3, 16- 
18,63,70,78]). Our goal here is to obtain some elementary formulas for the averaged Hamiltonian, which 
are probably new, and to give a global interpretation of the averaged motion basing on the geometrical and 
topological approaches to integrable Hamiltonian systems developed in [19,20,43]. We are also going to 
show that general result [71] gives probably the most complete statement about the averaging for H; it 
seems that the variables {Ii,ipi,y) are most convenient for the analysis involved 
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(3.4) 



To simplify further formulas, let us introduce the averaged potential v. Expand v into the Fourier series: 

v{xi,X2) = v{Q + yi,P + y2) 

= ?;(-v/27i sinc^i + j/i, •y/27i cos (/^i + 2/2) 

27ra2i 



^ Vk exp 

jk=(A:i,A;2)eZ2 

27r 



z(^/i:i(-\/2i7sin(/Ji + yi 



a22 



^2/1 COS(/5i + ^2)) 



+ ^2 



022 



2/1 COS (^1 + y2 



Now let us average the potential v with respect to the angle variable ipi : 

r2TT 



(3.5) 



1 /-^ 

v{Ii,yi,y2) = ^ 



V dip 



1- 



Taking into account expansion (13.41) and using Bessel's integral representation for the Bessel func- 
tions [52, no. 7.3.1], one can rewrite (13.51) as 



(3.6) 



v{h,yi,y2)= Yl vM^2h{kl + {2TiY{k2-kia2iY/al^) 
fc=(/ci,fc2)ez2 



X exp 



ikAyi 



27ra2i 

«22 



y2 + ik2 



27r?/2 

«22 



where Jq is the Bessel function of order zero. Using the spectral theorem, (13. 6t can be rewritten in a more 
elegant form: 

v{h,yi,y2) = Jo(\/-2/iA)v(yi,y2)- 
Here the operator Jo(\/— 2/i A) is a pseudo-differential operator [82]. Note that v is analytical with respect 
to Ji, because Jq is an even function. 

Let us formulate now our main result on averaging. 

Theorem 3.1. For any k > there exist Eq > 0, positive constants C and G, and a canonical change of 
variables 

P = T + eC/i(T, 0,^1,^2,6), Q = Q + e[/2(a',Q,yi,y2,e), 

2/1 = yi + eTyi(y,Q,yi,y2,e), y2 = y2 + eH^2(y,Q,yi,y2,e), 

defined in the domain Ji < k, e < eo (here and later Ji = ^(J*^ + O!^)), such that 

= ?{(Ji,y,e) + e~^/^g(y,Q,y,e). 
Here C/1,2. W^i,2. S are real analytic functions ofT, Q, y, and 

\Ui,2\,\Wi,2\,\S\ + \Vy5\<G, 

'K is a real analytic function ofJi and y. The functions Ui^2, W^i,2. S. ^ are periodic relative ^ with 
periods (ai, 02). In addition, we have the estimate 

:>{(Ji, y, e) = ^(Ji, y, e) + e''g{'y, q, e), ii-pi, y, e) = Ji + et;(Ji, y), 

where \g\ + |Vy(7| < M for some positive constant M independent ofe. 

Proof of the theorem follows immediately from the general result of Neishtadt [71] in the domain Ji > 
kq > 0. To include in our consideration the neighborhood of /i = 0, we need some its modification based 
on some special choice of generating function of the requested transformation. On the first step, one has to 
find a canonical change of variables (P, Q, y) ^ {P' , Q' ,y') that reduces the Hamiltonian to the form 

(3.7) H = H'{I[,y', e) + e^g{P', Q' , y' , e), 

where I[ = ((P')^) + (Q')^) /2> and g = 0(1) as e tends to 0. Let us try to find this change of variables 
using generating function S{P', Q, y'i,y2) = P'Q + y'iy2 + £s{P', Q, y'l, 1/2) from the equations 



(3.8) 



yi=yi+£ 



ds 
dy2 



y2 = y2 + e 



ds 
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Substituting (I3.8t into (I3.7t . one obtains the following condition on s: 

(3-9) QQpj-P'QQ=viP',Q,y'i,y2), 

where 

Q, y'l, = V(^^ {{P'f + , _ ^(Q + y'^^p' + ^2). 

Introducing polar coordinates /, if) by the equalities P' = \/27 cos ip, Q = \/27sin V', one can rewrite (I3.9t 
in the form = v. General solution of this equation can be written as s = fv dijj, but the function s can 
be non-analytical relative P' and Q; to avoid this, one should choose the integration constant in a special 
way, for example, 

s{P\Q,y'i,y2) = vi'^cos(p,V2I sin ip,y[,y2) dip 

+ J v{V2I cos (p,V2I sin ip,y[,y 2) dip 



p'=v^cosV, 
Q=\f2i sinV 

This procedure is then repeated, and the Neishtadt estimations [71] are used, see [24] or [45] for details. 
Note that on the first step described above one has H' = I[+ ev{I[,y'). 

We illustrate the above consideration in the special case of example (11.71) . Then we find 



(3.10) v{Ii,y) = AJo(V2/i)cosyi + SJo(/3V2/i)cos(/3y2). 

Properties of the Bessel functions [52, no. 7.4] give the estimations 

v{h,y) = A{1 - hi) cos yi + B{1 - ]^0^h) cos(/3y2) + 0{ll), as h ^ +0, 



v{h,y) = aJ — ^cos(y27i -7r/4) cosyi 

+ cos(/3y27i - tt/A) cos(/5y2) + 0(A~^^^)> as h ^ 00. 

VP 

4. Classification of the averaged motions 

4.1. A one-dimensional Hamiltonian system for the drift. Since the function M is a periodic function 
of y, it can be viewed as defining a Hamiltonian system in two different phase spaces, namely: 

(1) in Euclidean phase space $ = Rj, g ^ = ^ and 

(2) in the phase space <I) = M| q x T^. 

Obviously, these systems are integrable and equivalent to the equations 
(4.1) ai= const >0, 



(4.2) y = JVy:K(Ji,y,e), J 



-1 

1 



Eq. ( 14. U defines a family of "cyclotron" circles S'c(Ji), 3i £ [0,cxd), in the coordinates (IP, Q). 
The boundary Ji = of this family is the rest point T = Q = 0. For each fixed Ji, (I4.2t is a one- 
dimensional Hamiltonian system. The trajectories of (I4.2t are the connected components of the level sets 
of "K. Clearly, the solutions depend also on e, the action Ji, and other parameters, but we omit this depen- 
dence to simplify the notation. The system (I4.2t describes the slow drift of the centers of the "cyclotron" 
circles on the plane M^. 

It is now convenient to describe the trajectories using the topological and geometric theory of integrable 
systems developed in [19]. One may consider as a Morse-Bott function on three-dimensional surface 
J£(J', Q,y) = E [20, §1.8], or one may consider "K, for each fixed Ji, as a function of variables y. In 
the latter case, we suppose that Ji has only a finite number of non-degenerate critical points in the elementary 
cell (for generic potential v this property holds for almost all Ji), i.e. 'K is a Morse function on the torus 
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Figure 4.1. Drift in the x-plane 

covered by the plane M.y, and (I4.2t is a Hamiltonian system on the torus (or on its covering M^). We 
prefer this second point of view and find immediately a complete classification of its trajectories. 

Proposition 4.1. For any trajectory, 7, of (14. 2t on T^, one of the following assertions holds: 

(a) 7 is a closed contractible curve', 

(b) 7 w a closed non-contractible curve', 

(c) 7 an extremum point of !K; 

(d) 7 a saddle point of 'K or a separatrix. 

In the cases (a) and (b), we have periodic trajectories on the torus. This means, that for any trajectory 
7 = y (r) there is a T > (the period) and d = (^1,^2) G such that '^{t + T) = '^(t) + d ■ a, where 
d ■ a = diai + ^202. In the case (b) the vector d is non-zero; moreover, if both components di and d2 
of the vector d are non-zero, then they are relatively prime. In the case (a) we have d = 0. Obviously, 
the vector d is unique for each trajectory and does not depend on the choice of the trajectory on My covering 
7. Moreover, as different trajectories y(r) cannot intersect on the plane Ey, for fixed Ji exactly two non- 
zero vectors d with mutually opposite directions are possible; we fix one of them and denote it by (i(Ji); if 
necessary, we write ibd(Ji). The vector d defines the "average" or "main" direction of the motion {drift) on 
the plane (^1,^2) and it is one of topological characteristics of H', the meaning of this vector becomes even 
more clear if one considers the projection of the corresponding trajectory in the original space ^. onto the 
x-plane, see Fig. 14. II The ratio di/d2 is called the rotation number, see e.g. [§1.6] [20]. 

Definition 4.1. We call d{'ii) the drift vector^ of the motion. 

4.2. Tiie Reeb graph and the classification of the drift motion in non-degenerate case. Recall [20, 
Chapter 2] that it is possible to classify Morse functions on the torus by corresponding foliations, given by 
level curves, such that one has a foliation with singularities; the singularities are caused by critical points 
of the Morse function. There exist infinitely many topologically different types of such foliations which 
may be classified by their Reeb graphs. The complete theory of this classification is elegant but not trivial 
(see [19]), and we restrict attention here to the simplest situation assuming that 3i is a minimal Morse 
function on the torus T^, i.e. that "K has exactly one maximum point ^max and one minimum point ^niin 
(and hence two saddle points ^±). We put gmax.min := ^(^i, ymax.min), g± ■= ^(^i, ^i) and suppose that 
> g-. The classical motion for a fixed 3i is possible when g^in < < gmax- Recall that a minimal 
Morse function is called simple if 5+ > g- and complex otherwise. 

Thus assume first that for some Ji the function J{ is a simple Morse function. It is instructive to imag- 
ine that IK is a height function on torus (as shown in Fig. I4.2t : then the trajectories are the curves of 
constant height levels (compare with subsection 12. 4t . Consider the three intervals (^min, 5-)> {g-,g+), and 
(5+;5max)- If 5 £ (^mitu 5- ) U {g+,gmay.), then the set "K = g includes only one connected component, 
which is a closed contractible curve on T^, diffeomorphic to a circle. Hence we are in the case (a) of Propo- 
sition and this regime of motion is described by edges ii and of the Reeb graph, respectively: each 
point (denote it by gi or g^) on these edges corresponds to a contractible trajectory S'(fl'i,4, Ji) (see Fig. l4.2h ). 
The corresponding rotation number is equal to 0/0 and the drift vectors are d{ii^4) := (0,0). Each curve 
S{gi,4, 'Ji) C induces a set of closed trajectories Si{gi^4, Ji), I = {li, I2) G Z^, on the covering (plane) 
M|. We will discuss the numbering of Si{gi^i,3i) by / a little later. 

A closely connected notion appears in a more complicated situation in [75] 
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a) Reeb graph b) Level curves c) Action variables 



Figure 4.2. Characteristics of Morse function 



If 5 £ (9-^9+) thcri the set £ T'^ : "K = g} consists of two connected components, each of them 
being again a trajectory on the torus, diffeomorphic to a circle, but now they are both non-contractible such 
that we obtain a non-trivial rotation number di/d2. We in the case (b) of Proposition 14. II and each trajectory 
is characterized by the points g2 G 12 and 53 G ^3 on the two edges 12 and ^3 of the Reeb graph; denote these 
trajectories on the torus by S{g2.'i, Ji). They induce a set of open trajectories 5*^(52,3, Ji)> A; G Z, on the 
covering space M^. The numbering of these trajectories we also discuss later. The drift vectors d = d{i2) 
and d = ^(^3) corresponding to the edges 12 and ^3 have opposite directions; we use the notation d = 
for 12 and — for ^3. 

Points on the Reeb graph corresponding to the extreme levels g^ia and ^(max are called end points; they 
correspond to the stable rest points of the Hamiltonian system (I4.2t . The points corresponding to the critical 
levels g± correspond to separatrices, including the saddle points Thus each interior point of any edge 
defines a closed trajectory or a closed oriented curve on the torus where the orientation is given by the 
natural parameter t (the time). We may parameterize the edges of the Reeb graph by the variable g given by 
the value of IK; this defines one-to-one maps g : [gmm,9-) h, 9 '■ [9-, 9+] ^2, 9 '■ ^3> and 

9 ■ (5+,5max] H- 



4.3. Action variables and parameterization of the drift trajectories. One may also parameterize points 
on the edges of the Reeb graph by action variables 

1 r+^^ r ^.m f ^ M^){d{ij)-ah {d{ij) ■ a)i{d{ij) ■ a)2 



(4.3) h = ^ Mr)dMr 



2tt Att 



with sign prescribed by the natural orientation, where a is an arbitrary real number. 

As this definition of the action variable in the phase space is somewhat different from the one in M?, 
let us explain formula (14.31) . 

The second and the third term are present only for the edges i2 and ^3. The geometric interpretation of | J2I 
for zi 4 is standard: 27r|J2| is the square of the domain in My bounded by the corresponding closed trajectory 
(see Fig. 14.2b ). Of course, here the action variables do not depend on parallel transport of the coordinate 
system on My and on the choice of the closed curve on My. It is not difficult to show that J2 is positive for 
ii and negative for i^, that | J2I < ^22 = aiia22/(27r), and that J2 = in the end points (7min,raax of zi^4. 

The geometric interpretation of J2 for the edges Z2,3 is as follows. Denote by the straight line passing 
through the origin on the plane My in the direction of c? • a. Consider one of the lifts Sq = = ^(t)) on 
My of the trajectory S on Ty. Let us fix two points ^(cr) and ^(ct) + (^(^2,3) • a on this curve and project 
them onto L^, such that we obtain some curved trapezium. The square of this trapezium is equal to |27rJ2|, 
where J2 is defined by ( I4.3t (see Fig. 14.2b ). This interpretation allows us to derive some simple properties 
of the action variable. 

In contrast to the previous case we see that now J2 depends on the choice of a lift on the plane My of the 
trajectory (but it does not depend on this choice modulo aiia22/27r), and it also depends on translations of 
the coordinates. 

Let us fix next a certain continuous family of trajectories on the plane My corresponding to the edges ij. 
Using the geometric interpretation of 3 2 it is easy to show that 3 2 increases along each edge ij, such that we 
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obtain a parameterization of all trajectories on the torus by means of the action variable J2, associated with 
the Reeb graph. If necessary, we write to indicate that this action variable is associated with the edge ij . 
Obviously, J2 admits upper and lower limits along each edge. These limits depend also on 3i and given 

by 

j^+(Ji)= lim a^, a^(Ji)= lim a^, af/=^^(ai)= lim jf. 

g^g-~0 9^g++0 9^S±=F0 

One may calculate all the quantities as 

where 7=*= is a separatrix connecting the corresponding saddle point y ± ( J 1 ) with the saddle point y ( Ji ) + d 
on the plane M^. The numbers d^^^Ji) and J2 (^i) are well defined, whereas the numbers 32^^^^{'Ji) are 
defined only up to aiia22/2vr. If one fixes one of them, say, 3 2" then the others can be uniquely fixed 
by the "Kirchhof law" 

J^+pi) = + J^(Ji), 3l-^{3i) + dl-^i0i) = OtiOi) + ^aiia22. 

We fix ^2(^3,4) in this way; the choice of Jg" (^i) will be explained later. But for any choice of the action J2 
the following inequalities and equalities are true: 

(4.4) < 32 - ^r^^" (Ji) < ^aua22, 

Zn 

(4.5) + - + (^j3+p^) _ j3-p^)^ ^ j4-(j^) ^ ^«ll«22. 

Now we describe the numbering of the closed trajectories -^^(gi 4, Ji). We define the multiindex / as 
follows: Let us fix some extreme point ymin,max of !K in M^; we give the number / = (0, 0) to this point. It is 
clear that this choice determines the numbering of other extreme points by yf,,^ „• := yminmax + I ■ a, and 
the numbering of the corresponding trajectories Si{gi^4), depending continuously on J2 (they also depend on 
Ji, see subsection l4.5t . Indeed, if S'o(fl'i,4) is defined by the equation y = ^ (r, J), then the other trajectories 
are 5i (51,4, Ji) := y = ^ (r, J) + / • a. ' 

In contrast to the case (a), we enumerate the curves Sk by a single index A; e Z. Fix a vector / = 
(/i)/2) £ conjugate to d, i. e. difi + ^2/2 = 1, which always exists. Then we fix some open 
trajectory Sq corresponding to a certain point from the edge 12 and give it the index A; = 0. According to the 
"Kirchhof law" for the action variables, we have to give this index also to the full family of open trajectories 
associated with both edges 12, 13 depending continuously on the corresponding action variable 32- If these 
trajectories are given by the equation 6*0 ((72,3, ^i) : ^ = ^ (t, J), then the open ti"ajectories with index k are 
Sk{g2,3,:Si):'S=)^{r,3)-kJf-a. 

Thus we see that the trajectories of the system on the torus and on the plane are parameterized by action 
variables 3i and 32, indices Z G or G Z, and the edges of the Reeb graph; we include all these 
parameters to the notation in the next subsection. Finally, we have Figs. 14.2b and 14. 2b for the trajectories on 
My (generally speaking, the curves J2 aid J2 in Fig. 14.2b may coincide). 

4.4. The Reeb graph in degenerate cases. Now consider now the case when 'K{3i, •) is not a simple 
Morse function. 

There are two possible cases. Firstly, !K(Ji, •) can be a complex Morse function. Denote the correspond- 
ing value of 3 1 by 3\, and if necessary add the subindex a for numbering of these critical values. 

The regime of regular motion consists of contractible curves only, and the Reeb graph has the form 
described in Fig. 14. 3h . This graph may be considered as a limit of the previous case as — > 3+ such that 
and the edges i2 and ^3 contract to a common point. The action variables are sketched in Fig. 14.3b . and the 
phase picture for the trajectories on the plane is shown in Fig. 14.3b . 

Another limit case ('Fig. l4.4t occurs if g- — > ^min and 5+ — > gmax', all contractible trajectories disappear, 
and we are left with non-contractible trajectories only. In this case, "K is not a Morse function, but we can 
still assign a Reeb graph to this situation (see Fig. l4.4h ): we keep only the edges 12 and ^3, and the action 
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Figure 4.3. Characteristics of Morse function: The degenerated case I. 



variables and j| show a similar behavior (Fig. 14.4b ). We denote the corresponding values of the actions 
3i by and add, if necessary, the subindex a for numbering of these points. 

4.5. The description of the averaged motion in 4-D phase space. Using the above considerations we 
now represent the global structure of the classical motion defined by 'K under assumption that the domain of 
the motion on the half-plane (3i > 0, E) is separated into such (connected) subdomains, that the behavior 
of trajectories of the corresponding Hamiltonian systems for each of these subdomains is topologically 
equivalent and have the same rotation number (see Fig. Il.lt . As before, we call these domains regimes. 
The interior points of each regime correspond to closed trajectories on tori T^; the boundary of the regimes 
is formed by the critical manifolds of the function Ji and by the left boundary 3i = 0. 

A regime is called of boundary type and is denoted by Mr if it a certain part (of non-zero lenfth) of its 
boundary consists of extreme points of "K). On the plane My the corresponding level curves of the function 
Ji are families of closed trajectories, their preimages on the torus are contractible trajectories with ro- 
tation number 0/0, and they have no special direction. If we return back to the original four-dimensional 
phase space M^^, for each interior point in these regimes we get a family of invariant Lagrangian manifolds 
of Hamiltonian Ji. They are topological products of the cyclotron circles 5c (Ji) and the closed curves 
31(91,4,31), and they are diffeomorphic to two-dimensional tori (we call them Liouville tori). We parame- 
terize each point in by the action variables Ji and ^2, which belong to a certain domain on the plane M^; 
to simplify the notation we denote these domains also by M^. 

Thus each interior point J = (Ji,J2) G 'Mr indicates a discrete family of invariant Lagrangian tori 
A[(Ji,J2)^ in the original phase space M^^. numbered by a multiindex / = (/i,/2) G given by the 
numbering of the curves S'i(5i,4, Ji). 

If Aq (J) is defined by the equations 

(4.6) Xi,2 = ^1,2(3, ^), Pl,2 = A,2(J, f), 

where 99 = [ipi, ip2) are angle variables conjugate to J = (Ji, J2) S M^, then 

A[(J) = {x = X(J,v?) +/-a, pi = Pi{J,ip) + il-a)u P2 = P2i3 , if)} . 



'Of course, these manifolds depend also on e; we will include this fact into notation later. 
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The action variables do not depend on I, and they are defined by (I4.6t with d = 0. 

The remaining regimes are called interior regimes. We denote them by and use the symbol also 
for the domain of the corresponding action variables on the (Ji, J2)-plane and (Ji, E)-p[me. On the plane 
M.y, the level sets of are families of open curves with the main vector = (di, ^2), and their preimages 
on the torus T2 are non-contractible closed trajectories with rotation number di/d2- In the original phase 
space Mp^, these trajectories are covered by a discrete set of the families of invariant two-dimensional 

Lagrangian manifolds AJ^(J), 3 = (Ji, J2) £ ^r, of Ji. They are products of the "cyclotron" circles ScCJi) 
and the open curves S'fc((7i,4, Ji) and are diffeomorphic to two-dimensional cylinders; for brevity we call 
(J) Liouville cylinders. 

The cylinders A^(J) depend smoothly on J = (Ji, J2) S M^, their numbering by the index k is induced 
by the numbering of the curves Sk{gi,A,3i). Hence we have 



where the vector functions P and X define the Liouville cylinder Aq(J) by (I4.6t . 99 = tp2) are the angle 
variables, and J G M^. 

We can give formulas for the action variables Ji and J2 directly on the tori and cylinders by: 



where (p, x) belongs to A{,(J) or to A|^(/), respectively. Note that in the latter case J2 depends on index k, 
but in what follows we fix J2 by setting A; = in the definition and using this action for parameterization 
of all A^(J). We also assume that the families A[(J) and A^(J) depend smoothly on J = (Ji, J2) in the 
whole regime. Then, by fixing a certain cylinder from the interior regime associated with the edge 12 of 
the Reeb graph and giving it the index = 0, we determine by the "Kirchhof law" the action the cylinders 
corresponding to the edge ^3. 

It is well known that Lagrangian manifolds have the integer- valued homotopic invariants which are called 
Maslov indices and connected with the cycles on these manifolds. Obviously the Betty number (the rank of 
the cohomology group, or the number of basis cycles) of any Liouville torus A[ is equal to two, hence A[ 
have two Maslov indices and the Betty number of any Liouville cylinder A^ is equal to one, hence A^ has 
one Maslov index. Standard calculations lead to the following simple fact. 

Proposition 4.2. The Maslov index of the cycles 71.2 = (y'2,1 = const) on any torus A[(J) is equal to 2 
mod 4. The Maslov index of the cycle 7 = ((/?2 = const) on any cylinder A^(J) is also equal to 2 mod 4. 

Each non-degenerate point J = (Ji , 0) of the extreme boundaries of the regimes M,. defines in a 
degenerate torus, namely a closed trajectory, which is an isotropic manifold. In this case we have only 
cyclotron motion the drift is absent. The other non-degenerate boundaries between Mr and JAr define 
separatrices of a one-dimensional Hamiltonian system with the Hamiltonian 'K, which generates in ^. a 
two-dimensional invariant singular manifold of the Hamiltonian IK. 

The critical points on the boundaries of the regimes induce degenerate singular invariant manifolds; these 
manifolds together with their neighborhoods are called atoms. Generally speaking, there exist infinitely 
many topological types of atoms, one can find some classifications some of them in [19, 20]. We restrict 

our consideration to the simplest situation described in the previous sections, then all critical values on the 

1/2 

axis Ji are of the form . The Morse function "K changes its type, when Ji crosses these values (in our 
simplest case only rotation number changes; in more complicated problems a new Reeb graph can arise). 

The left boundary Ji = plays a special role in quantum applications, as it corresponds to the so-called 
low Landau bands. In the original phase space the corresponding trajectories belong the two-dimensional 
invariant subspace 3i = 0. This subspace presents only slow drift and the "cyclotron" motion is absent. All 
previous considerations concerning K remain vahd, but now the "limit" Liouville tori A[(0, J2) in ^p,x 



Al{J) = {x = X{I,ip)-kJf-a, p^ = Pi{3,ip) + kJf-a, p2 = P2{3,ip,e)}, 



(4.7) 
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just closed curves, and the "limit" Liouville cylinders A5J(0, J2) are open curves with drift vector d; these 
curves are isotropic manifolds also. 

The critical points (Ji = 0, g = g±) can be considered as zero-dimensional singular manifolds. In the 
degenerate case g+{0) = g-.{0), there exist only two boundary regimes; this case is not generic, but it 
appears in connection with the Harper equation (see below). The other degenerate case is gmm{0) = g^ (0) 
and 5max(0) = 5+(0). It appears, for instance, when v depends only on one variable. It seems that in this 
case one can separate the variables in the original spectral problem. 

The angle points J = (0, 0) in the left boundary correspond to the stable rest points of both the averaged 
and the original Hamiltonian "K and H. 

At last we remark that there is no reasonable definition of the Maslov index for an individual isotropic 
manifold A G M^n, if dim A < n, see e. g. [9, 13, 34, 68]. But if this manifold arises as the limit of a 
family of Lagrangian manifolds, one can associate with this manifold a Maslov index and make its use in 
the semiclassical approximation. Obviously, this applies to the problem under consideration. 

4.6. Example. We illustrate the considerations of this section by the example (11.71) - (I3.10t . Let us consider 
first H instead of "K, then 

5min = -(A\Jo{y^i)\+B\Jo{(iy^i)\), g,,,, = A\Jo{y^,)\+B\Joi(3y^i)\, 



g± = ± 



A\Jo{y^)\ - B\Jo{f3^/2D 



The first series of critical points, ^, a = 1, 2, . . . , is obtained from the equations 

Jo{3lJ = and Jo(/3J1,„) = 0. 
The second series, df ^, a = 1,2, . . . , consists of the solutions of the equations 

A\Jo{3lJ\ = B\M(33lJ\. 
On the plane (Ji, E), the boundary regimes are the sets (see Fig. ll.U 

Jl + e^min (Jl ) < ^ < Jl + £5- ( Jl ) , ^ia <^ia+l, 

and 

Ji + eg+CJi) <E <3i + e5max(Ji), < ^1 < ^L+i; 
the interior regimes are the sets 

Ji+eg-iJi) <E <3i + eg+ ( Ji ) , Oi < • 

Note that the rotation number changes when 3i crosses these critical points. 
The drift vectors of the interior regimes are 



1,0) iiA\JoU2J\j\> B\Jo{(3J2J\ 



(0,1) ifA\JoU2J\ )\<B\MP J 2Jl^ 



A simple calculation gives 



3l+{J,) = -JtiJi) = 23l+{J,) = 23l+{3,) = \ r'^^-Aog d^. 



(4.8) 3l^{3,)=3l^{3,) = -^--3\+(:3,) 



where 



and the sign in the exponential is such that T < 1. 

Using 'K gives a discrepancy O(e^) in all above estimations. 
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5. Almost invariant manifolds of the original Hamiltonian 

Definition 5.1. Let A = {p = P{ip,e), x = X{ip,e)} C Mp^. be either a two-dimensional Lagrangian 
manifold, diffeomorphic to a two-dimensional torus (or cylinder), or a smooth closed (or open) curve. Let 
C > 0. We say that A is an almost invariant manifolds of the Hamiltonian H{p, x, e) up to 0{e~'~^^^) if 

(5.1) i/|A = const + 0(e"'^/^), 

and if a vector u) exists, such that 

C5 2) {-^i'P + ^)7-^if + f^^' ^)) satisfies the Hamiltonian system up to 0(6"*^/^), uniformly 
in t e M. 

We call a family of two-dimensional almost invariant tori (respectively, cylinders) depending smoothly on 
action variables J G M almost Liouville tori (respectively, cylinders) of the Hamiltonian H. 

According to this definition, the manifolds A[(J) and A^(J) constructed in the previous section are almost 
Liouville tori or cylinders of the Hamiltonian H. 

Of course, not all the applications need a construction of this accuracy. In some cases, it is reasonable to 
construct an averaged Hamiltonian up to 0(e^+^) and to obtain "almost invariant" manifolds of H with 
a discrepancy 0(e^+^) (this means that in (I5.lt and (15.21) one has e^+^ instead of e"*^/^). In particular, H 
is the averaged Hamiltonian up to O(e^), and its invariant manifolds are almost invariant manifolds of H 
but only up to 0{e^). 

Let us fix a sufficiently small 5 and consider a boundary or an inner regime without ^-neighborhood of 
the boundary formed by separatrices. If the regime includes the left boundary Ji = 0, then a neighborhood 
of this boundary (with a neighborhood of the singular points removed) also belongs to the non-singular 
part of this regime. Let us denote these non-singular parts of Mr or by M.r,s or 3Vtr,(5> respectively. 
The corresponding regimes for to we denote by M,^^ and M^; they coincide with Jiir,s and Jiir,s up 
to O(e^), and they have the same drift vectors. Moreover, the corresponding almost invariant tori and 
cylinders, A['^ and A['^ have the same structure as A['^^ and A['^, and, in particular, their Maslov indices 
coincide. 



6. Semiclassical spectral series 

Now we are going to use the almost invariant manifolds described in the previous sections for constructing 
the spectral asymptotics for H, like it was done in subsection 2.5. Let us emphasize again that in even in 
the simplest multidimensional problems there is usually no global asymptotic formula for the spectrum; 
it is useful to divide the spectrum into several parts, such that the asymptotic behavior of the spectrum is 
preserved in each of these parts. These parts together with the corresponding formulas for the spectrum are 
usually referred to as spectral series. According the fundamental correspondence principle of the quantum 
mechanics (connected with names of Bohr, Sommerfeld, Einstein, Ehrenfest, Brillouin, and others), the 
classification of spectral series is connected with the motions of the classical Hamiltonian system [27, 46, 
60, 66, 67]. As a mathematical expression of this principle we use the method of the canonical operator 
developed by Maslov. 

Now we are going to show how the correspondence pPrinciple appears in the problem under study. 

Definition 6.1. A pair E) with h, e) e C^{E?) and E{h, e) e M is called a quasimode of the op- 
erator H with error 0{h^ + e^) if for any compact set 17 C there are positive numbers A and B 
satisfying \\{H - ^)^'||L,(f^) < Ah^ + Be^. 

From now on we fix some positive integer numbers K and L; in the rest of the section we describe the 
construction of quasimodes for H with error 0{h^ + e^). 

6. 1. Quasimodes associated with the almost Liouville tori. Consider a certain boundary regime Mr^s and 
the corresponding family of the almost Liouville tori A[(J, e). As it was mentioned above, the asymptotics 
of the spectrum of H is defined by means of the canonical operator. 
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For each fixed ^ > let us choose a discrete subset of the values of the action variables Ji and 02, by the 
rules . 

(6.1) 3, = 3'^(h) :={^ + fi)h, 

(6.2) 3^ = 3-{h):={^ + u)h, 
here fi and u are integer numbers such that 

(6.3) eM,,5. 

The conditions (I6.lt and (16. 2t are nothing but the necessary and sufficient condition for constructing the 
canonical operator on the tori A[ ( (/i) , (^) ; ■ 

Proposition 6.1. For any {fj.,u) £ x Z satisfying (16. 3t ?/jere ra'^f quasimodes [ip^f {x, h, e) , ''^ (^i s)) 
of the operator H with error 0{h^ + e^); f/j^^e quasimodes can be given by the equalities 



x'^f = i + o{h) G c-(A[(a^(/i),a^(/i),e)), 
= j{'^(j^(/i),a^(/i),e) +o(/i2). 



A// the functions ip'^'l^ belong to L^(M^) f/iey can be expressed through each other as (cf. (2.36)): 

(6.4) <f (x,/i,e) = V^'o'(^ - ^ • a,/i,e)e-i'2™; 

these functions are asymptotically localized near the projections vr^ A[ [3i{h) , J2 (^) ) o/f/je corresponding 
tori onto the x-plane as h tends to 0, /. e. 

lim V^^^f = Qasxi 7r,A[(a^(/i),a^(/i),e) 
One a/i'o has the estimate dist (^Er''^{h, e), spec = 0{h^ + e^). 

Remark. The numbers E'r ''^ as well as the functions depend also on K and L; now we omit this 
dependence to simplify the notation, but sometimes we will write -E'^'^- instead of E^'^ to emphasize this 
dependence. 

The proof directly follows from the general properties of the canonical operator; we describe the scheme 
of the proof in the Appendix. 

Denote by T,^^j^{h,e) the union of all possible points Er'^{h,e). We call this set the semiclassical 
spectral series up to 0{h^ + e^) corresponding to the boundary regime Mr,^. 

6.2. Quasimodes associated with the almost Liouville cylinders. Let us consider now a certain interior 
regime M^.^^ and the corresponding family of the almost Liouville cylinders A[(Ji, e)- 

To construct the canonical operator on these cylinders, we quantize only the action variable Ji by the 
rule (I6.lt . and J2 remains free. The integer numbers in (I6.lt are such that 

(6.5) (J^(/i),a2) gm,,5. 

Proposition 6.2. For any J2) G ^+ x satisfying (16. 5t there exist quasimodes 
{'iprkix, 32, h, e),Er{32, h, e)) of H with error 0{h^ + e^); this quasimode is defined by 

J/l qv- -A* 

Vr,k - ■^A-(j5'(h),J2,£)^''.fc' 

rr,k = ^ + 0{h)eC^[Al{3tih),32,e) 
E^^i32,h,e) =:K{3'l{h),32,e) +0{h'^). 
The functions ip'^i^ belong to L^^^(i?^); they can be expressed through each other by the equality 

(6.6) V^;f;-(x,a2,/i,e) =V^^(x + A:(J/)-a,J2,/i,e)e^^^™, 
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and enjoy the property 

(6.7) Vi^^(x + d-a,32,h,e)= VS^f^Cx, ^2, h, £)e'^^^^^-^d-ahx,-{d-aMd-ah/2) ^ 

where d is the drift vector associated with the section Mr,5, and the vector / = (/i, /2) is dual to d, i.e. 
difi + c/2/2 = 1- Theftinctions tpj^j^ are asymptotically localized near the projections 7rj;A£(jj'(/i), 32,^) 
of the corresponding cylinders onto the plane x as h tends to 0. 

For the numbers Er we have the estimate dist{Er{32,h, £:), spec H) = 0{h^ + e^). 

In contrast to the boundary regimes, the points Er{32, h, e) form vertical intervals in the domains Mr,5 
on the plane {E, 3i). Denote by e) the union of all these intervals. We call this set the semiclassical 

up to 0{e^ + h^) series corresponding to the interior regime M^.^^. 

6.3. Semiclassical spectrum. The union of the semiclassical series E^^^(^, e) and e) correspond- 
ing to all regimes will be called the semiclassical up to 0{h^ + e^) spectrum of the operator H. We denote 
this set by T,K,L{h, s). An example of the structure of the semiclassical spectrum is shown in Fig. 1.3. 

The statements of this section assert only that for arbitrary K and L one can find points from the spectrum 
of H in 0{h^ + e^) -neighborhood of the set T,K,L{h,s). Now a natural question arises: what kind of 
relationship between the exact and semiclassical spectrum of the operator H exist? In particular: does 
the whole spectrum of H belong to a certain 0(/7/^ +e^)-neighborhood of Sa l? At the moment the answer 
is unknown, but we hope that it is positive. Our expectations are based, in particular, on the existence of such 
relationship between the semiclassical and the exact spectrum for the Sturm-Liouville problem (Section 2). 

From the other side, each point of the semiclassical spectrum "asymptotically" has infinite degree of 
degeneracy in the sense that one can construct infinitely many linear independent quasimodes with the same 
semiclassical energy; if e is small enough, then there are no isolated points in As the exact spectrum 

does not have such properties (at least for rational values of the flux 77), it is clear that the semiclassical 
and the exact spectrum do not coincide. Nevertheless, the semiclassical spectrum gives some information 
about the exact one. Before discussing this question, in Section 7 we consider the situation when 77 is a 
rational number and try to understand the meaning of the magneto-Bloch conditions (1.3) and (1.4) for the 
semiclassical analysis. 

At last let us note that to construct semiclassical spectrum up to 0(/i^ -|- £^) it is enough to know only 
the averaged Hamiltonian H^, see Section 5. 

6.4. Relationship between h and e, and the widths of the Landau bands. Up to now, we did not assume 
yet any relationship between the parameters h and e, but in concrete physical problems such relationship 
may appear, say, /i = e", At > 0. The energy levels Er'^{h, e) and £^(32, h, e) depend on the parameters 
h and e in a regular way, and this means that increasing of K and L gives the correction to the energy 
levels corresponding to smaller numbers K and L. We cannot say the same about the functions ip'^f and 

ip'^f,, because the ratio like e/h appears in the formulas for these functions. This fact does not allow to use 
the Raleigh-Schrodinger perturbation theory based on the small parameter e (as it was done e. g. in [63] for 
the case of fixed h). 

The formulas for the semiclassical spectrum describe the well-known broadening of the Landau levels 
E^ = 'Jiih) (these numbers are infinitely degenerated eigenvalues of H for e = 0) implied by the ap- 
pearance of the electric field. For each /j,, the union of all the numbers Er{h, s) and '^^(32, e, h) for all 
possible values of /x, u, and 32 will be called the iJ,th semiclassical Landau band and denoted by L^(^, e). 

If e (s'max {3 1 ) — gmin (^1 ) ) < ^' ^^'^^ semiclassical Landau bands do not intersect, and one can calculate 
their widths: 

(6.8) diamL^(/l,£) = £(5max(J0 - 5min(a^)) + 0{e^ + h^). 
For the example (1.7) we have 

diamL^{h,£) = 2e(A|Jo(Y^)| +5|Jo(/3y^)| +0(/i2 + e2)). 
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For small and large values 0^ one can use the estimates for the Bessel functions (see Section 3); in 
particular, for large 3^ we have diam Lf^{h, e) « (Jj*)"^/^. For the example (1.7) we have 

diam L^{h,e) « e(A(2 - J^) + 25/3(2 - J^*)) for small J^*, 



diamL^(/i,e) « e{Aj —{\ cos( J2J^ - 7r/4)| + 




+ ^|cos(/3j2J^-7r/4)|)) for large 

VP ^ 

Note also that numerical considerations [51] show that the flux-energy diagram for the operator H for 
each Landau band in the plane {E, -q) looks like a butterfly ("Hofstadter butterfly"). Our assumptions about 
the smallness of h and e mean that we consider the asymptotics of the spectrum corresponding to the upper 
part of the Hofstadter butterfly. 

In this paper, we do not consider the structure of the spectrum in the neighborhood of the singular bound- 
aries; the standard semiclassical approximation does not work there. It is clear that this asymptotics depends 
of the type of singularity, and in some cases is based of the parabolic cylinder functions. We will study this 
question in forthcoming papers. 

6.5. Quantum averaging and the Harper-type equations. Basing on the Correspondence Principle, one 
can expect that there exist some quantum analogies of the averaging procedure (or, more generally, of the 
canonical transformation) in the classical mechanics. The study of such correspondence is developed now 
in several directions. 

One can consider the Schrodinger equation as an infinite-dimensional Hamiltonian system and use for 
its solution different nontrivial generalizations of the classical averaging methods (see [11, 58, 85, 86]). 
This view on the Schrodinger equation with operator H was used in [11], and this study is connected with 
non-commutative analysis. 

Another interpretation of the quantum averaging is based on the construction of an operator corresponding 
to the canonical change of variables, and this approach exploit the idea that a canonical transformation in 
the classical mechanics implies (with some accuracy) a unitary transformation in the quantum mechanics. 
This procedure approximately reduces the original spectral problem to the set of low-dimensional ones [54] 
and requires approximate solutions of some quantization problems. Such an approach for the problem under 
consideration gives the results formulated above. We have obtained these results by direct applying of the 
averaging and semiclassical methods; now let us try now to find a correspondence between the spectral 
problem for H and one-dimensional spectral problems basing on this interpretation of quantum averaging. 

Roughly speaking, according to [54], to construct the first nontrivial approximation (with respect to 
parameters h and e) to the solutions of the equation H"^ = E"^ one has to solve approximately the equation 



2' 9Q2 

where all the operations are ordered according to the Weyl rule. The operator in the left-hand side commutes 
with the harmonic oscillator ^{—h'^d'^/dQ^ + Q^); this observation lets us write the solutions of i6.9l in 
the form ^{Q,'^2) = ^/^(Q)^^^(y2)> where ^'/^ is the ^th eigenfunctions of the harmonic oscillator, and 
equation i6.9i is reduced to the family of equations 



(6.10) :K{3l{h),-ihj^,]}2,e)w^{Q,y2) = Ew^{'^2), ^i G Z+. 



In particular, if we take H instead of "K for the example (1.7), then the equations ( 16.101 ) will have the form 

+ BJo{(5^J^^)cOS^5^2W^{^2) = ^^W^{^2). 
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i.e. they form a family of the Harper equations. Using this analogy, the equations (I6.10t are usually called 
the Harper-type equations. Approximate solutions of each of them can be found using the usual one- 
dimensional WKB method. The description of the spectral asymptotics for equations of such kind using 
the Reeb graphs was obtained in [39]. Thefore, from this point of view, each semiclassical Landau band 
is described by a certain Harper-type equation; this equation depends on the band, and, therefore, different 
Landau bands can have different asymptotic structure. As follows from the preceding, the asymptotics of 
//th Landau band can be obtained by the quantization of the Reeb graph for the function J{ y i , ^2, e) 
considered as a function on the torus / (oi, 02), what can be easily seen from Fig. 1.3. 



7. The spectral asymptotics in the case of rational flux 

Consider now the case when the flux rj := a22/h is a rational number, rj = N/M, where N and M are 
mutually prime integer numbers and M > 0. As mentioned in Section 1, in this case to each point from 
the spectrum of H one can assign a family of eigenfunctions satisfying the magneto-Bloch conditions (1.3) 
and (1.4). Crearly, the functions ip^,'^ and ^^'^ do not satisfy these conditions, vbut, like in subsection 2.5, 

we can use them as a base for construction of quasimodes j = 0, . . . , M — 1, satisfying (1.3) and (1.4) 
(cf. [73, 87]). For convenience, we call the set of M quasimodes satisfying the magneto-Bloch conditions 
as n family of magneto-Bloch quasimodes. We expect that this procedure will improve a detalization of the 
asymptotics to the spectrum. 



7.1. Magneto-Bloch quasimodes corresponding to the boundary regimes. Let us consider a certain 
boundary regime Mr,^ and the corresponding semiclassical spectral series Tl^ ^{h, e). 

Proposition 7.1. For any (gi, 92) G [0, ^/M) x [0, 1) and any (J^, Jg) £ ^r,s there exist exactly M families 
of magneto-Bloch quasimodes of the form {cf. (2.18)) 

= ^Cf'^■(gl,g2,M€o(^-^«>/i>£)e-*"^'^"^^ 
(7.1) iez2 

s,j = 0,...,M-l, 

such that all the functions ^^fl%u are linearly independent. The coefficients C^''' {q) in ( I7.lt may be chosen 
in the following form {cf. (2.38)): 



( 

exp 



(7.2) c;;f;^(gi,g2,/i) 



2Tri{qili + q2n) + 27rzr//ij - zr//2a2i/2 

ifl2+j-s + nM = 0, n G Z, 
0, otherwise. 



Here Z = (/i, I2) G Z^, the index s = 0, . . . , M — 1 indicates the family of magneto-Bloch quasimodes, the 
index j = 0, . . . , M — 1 indicates the number of a member in each of these families. 

To prove this Proposition one has to substitute the sum ( 17.11 ) into (1.3) and (1.4), to equate the coefficients 
of iIjq''^'^{x — I ■ a,h, e), and to study the infinite linear system obtained. 

We see that for each fixed value of the quasimomentum q there are magneto-Bloch quasimodes 
corresponding to the same spectral value; in the other words, we have a degeneracy of degree M^. We try 
to give an interpretation of this fact in Section 8. 

Let us describe the structure of the functions ^ll^u'' ■ If the number is small enough, then the asymptotic 
support of each of them consists of family of annulas. These annulas form strips separated by array of M — 1 
"empty" strips (see Fig. I7.lt . where the corresponding quasimodes have order 0{h°°). If -q tends to an 
irrational number, then M — > cxd, and only one strip is kept. This means, probably, that each of generalized 
eigenfunctions in the irrational flux case is asymptotically localized in such isolated strips. The diameters 
of the annulas depend on Jj' (i.e., on the index of the Landau band); if the number Jj* is large, then these 
annulas may intersect, and, probably, can cover the whole plane 
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Figure 7.1. Asymptotic supports of magneto-Bloch quasimodes "^^^f^u on the x-plane 



7.2. Magneto-Bloch quasimodes corresponding to the interior regimes. Now consider the quasimodes 
associated with the almost Liouville cylinders. Like in the previous subsection, the functions do not 
satisfy the magneto-Bloch conditions, and we again use them as a base for constructing magneto-Bloch 
quasimodes, i.e. we put 



(7.3) 



%{x,h,£,q) = ^Cl{qi,q2,h) i>'^{x,'J2, e, h), q = (^1,92)- 



To obtain the expression for the coefficients C^, let us substitute the expressions (17. 3t into (1.3) and (1.4), 
take into account the property ibJl . and then equate the coefficients of + k{Jf) ■ a, 32,£, h) for all 

k. This system has the following form: 

\ i 1 
Cl exp 27r-J2 + 2-Kiki] - -i{2-Kdi + 021^2)^2 



(7.4) 



(■j+d2) modM 



exp 



L'< 



d2+j 
M 



2TTidi{qi — jrf) — i7]d\a2i/2 + 2TTiL'q2 



<{-L' + l)-—, L'eZ, 



Cl_^ exp i{k - l)/ir/(-27r/2 + a2i/i 



C. 



(i+/i)modM 



exp 



L" < ll±i < {-L" + 1 



2TTif2{qi - jv) - 021/2 + 27riL"q2 
1 

M' 



L" G 



(here and later by A mod M we mean the remainder after the integer division of A by M). Obviously, the 
system depends crucially on the drift vector d of the regime. Let us consider first the case d = (±1,0), then 
the system for the coefficients takes the form 



(7.5) 



Ci'^e^2n^(3^+ka22) ^ (Jh±^-2^t{qi-jv)^ ^ ^ j = 0, . . . , M - 1, 

keZ, j = 0,...,M-l, aM_l = e2-'^ aj = l, j ^ M - 1, 



where the index it corresponds to d = (±1,0). We see that all the coefficients are uniquely determined by 
arbitrary chosen numbers Cq'^, j = 0, . . . , M — 1, and we have, therefore, at most M families of magneto- 



Bloch functions. For sth family, s = 0, . . . , M — 1, set Cq'-''^ 



It is easy to see that the equality 



(7.6) 



^2 = ^2 =^\J^^1^ 



n 



where is the action variable corresponding to the drift vector d = (±1, 0), is a necessary condition for 
the existence of solutions for M .51 . Obviously, the coefficients C^'-'' can be obtained from one set, say, 
C^'-''^, by the shift of the index j (and this means that really we have only one family of magneto-Bloch 
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Figure 7.2. The structure of the asymptotic support for the magneto-Bloch quasimodes 
(a) the drift vector is (1, 0); (b) the drift vector is (0, 1) 



quasimodes). Therefore, the resulting coefficients can be chosen in the form 

' ^tvk^a2,/2+2ntnq2 ^ if j k + uM = 0, 

n = n^N + hM, n G 
0, otherwise, 



(7.7) 



where N is an integer number such that for some M G Z one has NN + MM = 1. 

The expressions for the coefficients for the cases d 7^ (±1, 0) are rather complicated. At least, it is 
clear from the system (17 .4t that all the coefficients are non-zero in this case. From this point of view, the 
case d = (±1, 0) is the most useful one, because the asymptotic support of any magneto-Bloch quasimode is 
minimal in this case. From the other side, a situation with an arbitrary drift vector can be reduced to the case 
d = (±1, 0), if one turns the coordinates, applies a gauge transformation, and transform the magneto-Bloch 
conditions accordingly. But it is important to emphasize, that this reduction is not global, because the drift 
vector can jump when passing from regime to regime even in the simplest cases, and to obtain reasonable 
formulas for the magneto-Bloch quasimodes one has to apply different transformation for different regimes 
(or, respectively, these transformations depend on the Landau band), and there exist no "globally good" 
coordinates. 



8. Discussion and heuristic estimate of numbers of subbands 

8.1. Relationship between the true and the semiclassical spectra. Dispersion relations. Let us give a 
qualitative description of the semiclassical spectrum in the rational flux case basing on the considerations of 
the previous section. 

Let us consider a fixed semiclassical Landau band with index /i. Its part corresponding to a boundary 
section is discrete, and each point Er'^ is M^-fold degenerated. The part of the band corresponding to an 
interior regime is continuous, but now each point Er{32, h, e) is only 2M-fold degenerated as Er{32, h, e) 
can coincide for different edges of the Reeb graph. (Here we recall that we consider the situation with the 
simplest Reeb graphs, see §4.2; otherwise these estimations are estimations from below.) From the other 
side, if the number rj is rational, then the true spectrum of H has band structure and each point of the true 
spectrum is M-fold degenerated. Let us try to give an interpretation of the ambiguous degeneracy of the 
points Er'^ . 

The existence of the isolated points (J^*, Jg) can be explained then as follows. Our previous considera- 
tions have a non-avoidable error 0{h°° + e°°). Therefore, we can expect that these isolated points really 
approximate subbands of width 0{h°° + e°°). The presence of Af^-fold degeneracy of these points proba- 
bly means that in a neighborhood of each such point there are M true spectral subbands (minibands) of the 
operator H. (Additional arguments can be given basing on the following idea: the formulas (17. 2t realize 
a representation of the magnetic translation group on the space of asymptotic eigenf unctions; as all such 
representations are M-dimensional [95], small variation of parameters leads to the splitting of each energy 
level into M numbers.) Enumerate these bands by an index s = 0, . . . , M — 1, then the dispersion relations 
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Erjs{qi, q2, h, s) in all these minibands look as follows: 

Ei^;,-{qi,q2, h, e) = E>^.'J,,^{h, e) + Oih"^ + e^) 

for any positive K and L. Therefore, the approximation used does not give expressions for the dispersion 
relations, and the presence of the separation of bands cannot be found rigorously. Nevertheless, this splitting 
is supported by the analogy with the one-dimensional periodic problem (subsection 2.3). Assume that there 
exists a function lip!^^ such that the true magneto-Bloch functions can be represented as 

where the coefficients Cf'-^ are defined in Proposition 17.11 Similarly to (2.23) we obtain that if and ^'2 
are (generalized) eigenfunctions of H with eigenvalues Ei and E2, then 



(8.1) E1-E2 



3f? f 

dD 



"St '^i'^2dxidx2 

where A = (— X2,0) is the vector potential of the magnetic field, ds = {dxi,dx2), and I? C is any 
domain with boundary dD. Assume that the asymptotic support of tp^^ belongs to the unit cell generated 
by the vectors ai and 02. Let us choose this unit cell as the domain D. Put ^'^^2 = ^ll^u'^i^, q^^'^, h, e) and 
Ei/2 = Er'siq^^"^, h, e). Substituting all these expressions into (I8.lt . one obtains (cf. (2.25)): 

(8.2) Et'^tf,ft,E)-E^i'■(9^/.,E) 



h,n=0,±l 



where 



Pn,h 









ds 






'^dxidx2 



I = {h,nM). 

In the interior regimes, we have the following dependences of the energy on the quasimomenta (semi- 
classical dispersion relations): 

E^{q,h) = E^±{q,h,e) = Ei^{j^{q,n^,h),h,e). 

Consider the case with drift vector (±1,0), then these functions depends essentially only on qi, and the 
dependence on q2 is absent up to 0(/i^ + e^). As some of these functions increases in qi and others 
decreases, for some critical values of qi = ql one has 

E^+{ql,h,e)=E^-{ql,h,e) 

(see Fig. 18. U . For the example (1.7), these points correspond to the values ^2^^ = h{2n + 1)/{2M). We 
can expect that these points together with the "end" points J^/^ = hn/M are 0{h°° + -approximations 
of the gaps in the spectrum of H, see Fig. 18.21 This expectation is based, in particular, on the analogy with 
the one-dimensional periodic problem (subsection 2.5), cf. [39]. If the interior regime has the drift vector 
other than (±1, 0), then these semiclassical dispersion relations depend on a certain linear combination of 
qi and q2. 

Our hypotesis about the structure of the Landau bands are illustrated in Fig. 1.5. 

Note that in our spectral estimates we have used the almost invariant Liouville tori and cylinders. It is 
known that even an exponentially small correction, which is kept after applying the averaging procedure, 
can destroy some of these objects, and non-Kolmogorov sets may appear. These fact implies the following 
question: what do non-Kolmogorov sets mean for the exact spectrum of the operator H, in particular in 
situation when the flux 77 is rational? It seems that the answer cannot be given using the additive asymptotics. 




8.2. Heuristic estimate for the numbers of subbands. If the hypotesys of the previous subsection is true, 
then it is possible to count the number of spectral subbands corresponding to a fixed (semiclassical) Landau 
band. Let us consider a certain fixed value 3^. 

For the edge ii of the corresponding Reeb graph, the number of the quantization points J2 is equal 
approximately to — J^'^/h, and for the edge ^4 this number is equal to J2 Each of these points subbands, 
so the end edges ii and 14 give us approximately (modulo singular boundaries effects) 

n 011022 

bands. 

The expected numbers of subbands implied by the edges ^2/3 depends may depend on the symmetry prop- 
erties of the potential v. For the example (1.7), we obtain approximately (again modulo singular boundaries 
effects): 

rT2-f' rT2 — a2-|- n^i — 



2M ^ . ^ = 47riV 



h 011022 
bands. 

Therefore, the total number of the subbands is approximately equal to 

aH- 't4— _|_ q2+ q2— I q3+ q3- 

J2 ''2 ' ''2 ''2 ' ■'2 ■'2 



2t:N 



011022 
t2+ a2— a3+ a3- 



(we have used the symmetry property J2 ~ ^2 ~ ^2 ~ -^2 ) ^^'^ latter number is precisely equal to 
according to the Kirchhoff law (4.5) for the action variables. 

8.3. Correspondence to difference equations. In conclusion, let us discuss the relationship between the 
problem under consideration and the difference equation (subsection I6.5t in the rational flux case. It is 
known that the spectrum of a rational Harper-like operator with flux rj = N/M consists of N bands [39]; the 
presence of M -grouped bands corresponding to finite motion (boundary regimes) was studied numerically 
for some cases [39]. Therefore, we can expect, that our hypotesis about the spectrum of H, in particular, 
the estimate of the number of subbands in each Landau band is not connected with the simplicity of the 
potential. 

Appendix A. The canonical operator and spectral estimates 



As it was noted above, the asymptotics of the spectrum can be found using the canonical operator. Let 
us remind some basic properties of the canonical operator; a more detailed constructions can be found, for 
example, in [60] or [68]. 
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Let A be a closed Lagrangian manifolds without boundary in the space Mp" . The canonical operator %\ 
corresponding to this manifold maps from (7°° (A) to C°°(IR"), and for any function / G C°°(A) one has 
Kf = outside certain (5-vicinity of the projection tt^^ The canonical operator on A can be constructed iff 

1 / , / Ind7\ ^ 

— (Ppdx = n[n-\ — , n ^ £j, 

271 J ^ \ 4 / 

for each basis cycle on A, where Ind denotes the Maslov index. Applying this consideration to the tori A[ 
and the cylinders AJl we obtain the quantization condition (I6.lt and (I6.2t for the tori (they have two basis 
cycles), and (I6.lt for the cylinders (only one basis cycle). 

The following commutation formula is one of the most important properties of the canonical operator. 

Proposition A.l (Commutation formula). Assume that A is an invariant Lagrangian manifold of a Hamil- 
tonian system for a certain Hamiltonian H, and that on A there exists a volume form which is also invariant 
under the Hamiltonian system. There exists a sequence of differential operators {W}'^^, 

W : C°^(A) ^ C°°iA), 

with smooth coefficients such that for any function ip £ C°° (A) and any number G N one has 

N 

(A.l) H%K<p = %x[^{ihyR^v) + 0{h^+^). 

j=0 

In particular, is the operator of multiplication by the scalar function H\\, and 

^1 _ d ^ dH d dH d 
dt dp dx dx dp 

Here H is the Weyl quantization of the Hamiltonian H. 

Now let us try to construct an approximate solutions of the equation H"^ = up to 0{h^ + e^) basing 
on the commutation formula and almost Liouville tori A[. Assume that the conditions (I6.lt and (I6.2t are 
satisfied. Let us find requested solutions in the form 

L-l L-l 

^ = XAru, n=^(i/iynj, UjeC°°(A[), S=^(i/iyAj. 

j=0 j=0 

Applying the commutation formula, we obtain 

L-l 

{H-E)^ = {Y, E - ^^>4 + ^^^^ + 

n=0 s+j=n 

(the presence of the term 0{e^'^^^) is implied by the fact that the manifolds A[ are not invariant under H, 
but only almost invariant). We request that the expression in the curly brackets vanishes at least up to 0{e^) 
for some K > 0. 

For n = we obtain the equation (i? |a[ — Ao)tto = 0, and we can put Aq = ^Iaj" and uq = 1. 
For n = 1 we obtain {d/dt + Ai)no = and we set Ai = 0. 
The equations forn > 2 have the form 

d 

-^Wn-l = y^jR^'^ - \n-j)Uj 
j=0 

(they are called homological equations). Let us show that all these equations can be solved up to 0(e"^), 
where m is arbitrary positive number. To do this, let us note first, that the operator d/dt on each torus can 
be written as 

d a d dJ{ 

TtT^^I^ ^"^2^ , /2 = ;^;| ■ 

dt Oipi 0(p2 c)Ji/2 A[ 

It is important for us that = 1 + 0{e) and uj2 = 0{e); both these numbers do not depend on (pi and (p2. 
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Rewrite all the homological equations in a common form 
(A.2) ±f = E + g, gGC^m. 

Let us expand all the functions into their Fourier series: 

(fci,fe2)ez2 (fci,fc2)ez2 
Substituting these series into (IA.2t . we obtain /orma/Zy: 

^'^''^ = hu;,l hu;2^'^^''' ^ = -50,0. 
Note that these coefficients fki,k2 can be very large because of the denominator, and the function / is, 
generally speaking, not defined. To avoid this obstacle, let us use additional estimates. As g £ C°°, its 
Fourier coefficients decay very fast, and for any a > there exists positive numbers C(q) and N{a) such 
that |5fci,fe2l — ^('^)/(l^il + l^al)'^ as \ki\ + \k2\ > N{a). Let us put a = 2m. 
Introduce a set Q{e, a) as follows: 

Q{e,a) = [{ki,k2) G : |fci| + |A;2| < N{a)^ 

u[{ki,k2)ez^ : 1^21 <^}. 

Clearly, \kiLOi + k2i02\ > 1/2 as (fci, k2) G Q{s, a) and e is small enough. 
Now set 

G= E 9k,Me'^''''''^''^'\ g = g-G. 

(fci,fc2)eQ(£,a) 

The equation df /dt = G can be solved in Fourier series, and the function g gives a discrepancy O(e^), 
because 

~9k,M<C{a)/{\ki\ + \k2\r <C{a)e''. 
Therefore, we can construct a function u G L^(M^) and a number £' G M such that \\u\\i2 > c > and 
\\{h - E)u\\l2 = 0{h^ + e^) ash,e^ 0. Then 

(A.3) dist(E,spec^) < ll^^~f^'"ll =o{h^ + e^). 

\\u\\ 

The same procedure can be applied to each of the quantized cylinders A^, but as a result we obtain a 
function u G Lj^^ and a number E £R such that 

^x(x + d • a, /i, e) = n(x, /i, e)e^(2"^^-(*")^"i-^('^-")i('^-'^)^) , 

(A.4) 

lkllL2(n) > c> 0, and \\iH - E)u\\l2^u) = 0{h^ + e^), 

where 

Ii= {x = Ti{d-a) + T2J{d-a), riG[0,l], ra G M}. 
As such function u does not belong to L'f^^i^), one cannot apply the inequality (IA.3I) directly. 
Proposition A.2. Le? a function u and a number E satisfy the conditions \A.A\ . then dist(i?, spec -ff) = 

Proof. Introduce new coordinates 

a p\ f a\ d ■ a 



y = Ax, A 
and a function 



p aj' \I3) \d-a 



S{y) = li-apyf + a^yl + 2p^yiy2). 
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Define a unitary operator tj in L^(R^) by the rule 

f{x)^g{y)=e-i^^y^f{A-^y); 
it is easy to see tliat U is well-defined also on Lf^^{R^). Now we set H = U HU^^, i. e. 

1/ d \21/ 5\2 _-, 

As the operator U is unitary, the spectra of H and H coincide. 
Put if = Uu, then 



fivi + \d- a\,y2,h,e) = (p{yi,y2,h,e)e h 



Denote 



Note that 



n. = {(y 



i,y2) G 



-s\d ■ a\ < yi < s\d ■ a\> , s G Z, 



L2(n,) - Vs||/||i2(nj 



for any function satisfying f{yi + |d • o|) = e*°/(yi, 1/2), a G M; in particular, this holds for / = 99 and for 
Choose now a smooth function e(^) such that 

< e(e) < 1, 

e(.^) = 1 as ^ G {—\d ■ a\,\d ■ a\), 
e(^)=0 as^ ^ (-2|d-a|,2|d-a|), 

and choose a constant Cq such that 

|e| + |e'| + \e"\ < Cq. 

Put 65(7/1, 2/2) := e(yi/s). Now we have the following chain of equalities and inequahties: 

-y/sdist [E, spec H) \\ip\\ < dist [E, specH) ||es¥'||^2(n^') 

... 1 o . .... 

< 



1 de 

'H-E){esV) < es<^>- -h'^Aesip-h^{Ve\Vip) -ihx2-^ 



< \\es<^\\ + ^h'^\\Aes(p\\ + h'^\\{Ves\V(p)\\ +h 



oyi 



+ 



h'^Co^/su dtp 



+ 



X2ip 



L2(ni) 



s Il5y2"^'("^) 
Tending s to +00, we obtain the inequality 

dist(£;,specF)||v9||^2(n^) < Co\/2||(^ - E)ip\\^^^^^y 

Now one only has to use the conditions (IA.4I) . 



□ 



Note that in the problem under consideration there exist non-Lagrangian (but isotropic) invariant mani- 
folds of the averaged Hamiltonian. The construction of spectral series corresponding to such manifolds is 
studied in [23]. 
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